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Abstract 



micrOMEGAs is a code to compute dark matter observables in generic extensions 
^ of the standard model. This new version of micrOMEGAs is a major update which 

includes a generalization of the Boltzmann equations to accommodate models with 
asymmetric dark matter or with semi- annihilation and a first approach to a gener- 
alization of the thermodynamics of the Universe in the relic density computation. 
' Furthermore a switch to include virtual vector bosons in the final states in the 

annihilation cross sections or relic density computations is added. Effective opera- 
Qh' tors to describe loop-induced couplings of Higgses to two-photons or two-gluons are 

introduced and reduced couplings of the Higgs are provided allowing for a direct 
comparison with recent LHC results. A module that computes the signature of 
DM captured in celestial bodies in neutrino telescopes is also provided. Moreover 
^ , the direct detection module has been improved as concerns the implementation of 

I the strange "content" of the nucleon. New extensions of the standard model are 

' included in the distribution. 

p 

O ■ 1 Introduction 

cn _ 

Cosmological and astrophysical observations have established that the matter content of 
the Universe is dominated by dark matter (DM). However the nature of this dark matter 
^ . still needs to be elucidated, therefore searches for dark matter signals are pursued actively. 
c3 . Dark matter observables include the relic density, the elastic scattering cross section on 
nuclei, the indirect signatures of DM annihilation in the galaxy (leading to photons and 
antimatter) or DM capture in celestial bodies (leading to neutrinos). 

For all these observables, experiments are currently taking data and improved results 
are expected. The PLANCK satellite has recently determined the DM relic density with 
a precision of 2% [1]. Direct detection experiments such as DAMA [2], CDMS [3] Co- 
GeNT [4] and CRESST [5] have found events that could be compatible with DM in the 
range 5-30 GeV, these are however incompatible with the upper limits on spin independent 
DM interactions with nuclei set by XenonlOO [6]. Furthermore COUPP [7], Picasso [8] 
and XenonlOO [6] have set limits on spin dependent interactions. Larger direct detection 
projects with ton scale detectors are planned, e.g. Xenon, and should improve by more 
than one order of magnitude the current sensitivity. Indirect detection experiments are 
measuring the spectrum of cosmic rays in different channels. FermiLAT has measured 
the photon flux from the Galactic Center and from Dwarf Spheroidal Galaxies, the lat- 
ter have allowed to probe the canonical DM annihilation cross section for DM lighter 
than 40GeV [9]. Furthermore an anomaly corresponding to a gamma-ray line at 130 
GeV has been extracted from the FermiLAT data [10]. FermiLAT [11], PAMELA [12] 
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and AMS [13] have measured the positrons fluxes, showing an anomaly at high energies 
while PAMELA data on anti-protons are in perfect agreement with expectations from the 
background [14] . AMS will pursue the measurement of antimatter (positrons, antiprotons 
and antideuterons) and photon fluxes. Despite some anomalies, none of these searches 
lead to conclusive evidence for DM. Finally neutrino signals from annihilation of dark 
matter captured by the Sun are being probed by neutrino telescopes such as SuperK [15], 
ANTARES [16] or ICECUBE [17]. 

On a different front, the LHC has discovered the Higgs boson [18, 19] and has searched 
directly and indirectly - for example through their effect on B physics - for new particles 
that are predicted in extensions of the standard model with a DM candidate. Matching 
the properties of the Higgs particle imposes strong constraints on the parameter space of 
DM models. Furthermore the LHC is searching for DM signals in events where a pair 
of DM particles are produced in association with a jet or a photon [20, 21, 22], these 
searches provide the best limit on effective DM interactions with nucleons for light DM. 
Confronting the predictions from a variety of particle physics models which possess a 
stable weakly interacting massive particle (WIMP) with all these observations will help 
to corner the properties of a new cold dark matter candidate or show that one has to turn 
to other candidates that act rather as warm dark matter. 

The interpretation of the abundant recent and upcoming data from these different 
types of experiments requires tools to compute accurately the various signals for DM 
interaction and this in the context of different particle physics models. micrOMEGAs is 
a code designed to predict dark matter signatures in various extensions of the SM which 
possess a discrete symmetry that guarantees the stability of the DM particle.^ This code 
includes some of the common extensions of the SM : the minimal supersymmetric standard 
model (MSSM) and its extensions (NMSSM,CPVMSSM), the little Higgs model, the inert 
doublet model and incorporates various constraints from collider physics and precision 
measurements for some of these models. Furthermore the code can accommodate user 
defined models. In this new version we extend the functionalities of micrOMEGAs in three 
main directions: 

• generalizing the Boltzmann equations to include asymmetric dark matter and semi- 
annihilations. The DM asymmetry is taken into account when computing di- 
rect/indirect detection rates. 

• incorporating loop-induced decays of Higgs particles to two-photons and two-gluons, 
and computing the signal strength for Higgs production in various channels that can 
be compared to results from LHC searches. Limits on the Higgs sector can also be 
obtained through an interface to HiggsBounds 

• adding a module for neutrino signature from DM capture in the Sun and the Earth 

In addition the following improvements to the previous version are incorporated 

• Annihilation cross sections for some selected 3-body processes in addition to the 
2-body tree-level processes. The 3-body option can be included in the computation 
of the relic density and/or for annihilation of dark matter in the galaxy. 

^DarkSUSY [23], Isared [24] and Superlsorelic [25] arc other public codes that compute the signatures 
of dark matter. These codes are confined to supersymmetric models. 
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• Possibility of using different tables for the effective degrees of freedom {geff(T), heff{T)) 
in the early Universe 

• Annihilation cross sections for the loop induced processes 77 and 7^" in the NMSSM 
and the CPVMSSM 

• New function for incorporating DM clumps 

• New function to define the strange quark content of the nucleon 

• The LanHEP [26] source code for new models is included 

• New models with scalar DM are included (Inert doublet model [27] and model with 
Z3 symmetry [28]) 

• New implementation of the NMSSM which uses the Higgs self-couplings and the 
particle spectrum calculated in NMSSMTools_4.0 [29, 30] 

• New versions of spectrum generators used in the MSSM (Suspect_2.4.1) [31] and in 
the CPVMSSM (CPsuperH2.3) [32] 

• Extended routines for flavor physics in the MSSM 

• Update in interface tools to read files produced by other codes, this allows easy 
interface to other codes 

micrOMEGAs3.1 contains all the routines available in previous versions although the 
format used to call some routines has been changed. In particular global variables are 
used to specify the input parameters of various routines. All functions of micrOMEGAs are 
described in the manual, manual3.pdf, to be found in the man directory, see also the 
lecture notes [33]. 

In this paper we first describe in section 2 the generalization of the relic density 
routine to include asymmetric dark matter and semi-annihilation processes as well as 
more general thermodynamics of the Universe. Section 3 describes the new module to 
compute the neutrino-induced muon flux from DM captured in the Sun or the Earth. In 
section 4, we give a brief description of improvements on the computation of the direct 
and indirect DM signals, in particular to take into account possible DM clumps. In section 
5 we describe new facilities to compute more precisely the partial widths of the Higgs and 
incorporate LHC results on the Higgs. A brief description of SUSY models is provided 
in section 6. Section 7 contains a description of the functions that allow to work with 
external codes. Finally sample results are presented in Section 8. 

2 Generalization of the relic density computation 

The Boltzmann equation that describes the number density of DM that is used to solve for 
the abundance of DM today and obtain its relic density is based on several assumptions. 
First, it is assumed that the DM candidate is either a self-conjugate particle or that there 
is no DM anti-DM asymmetry, and that there exists a discrete symmetry such as R-parity 
in supersymmetry that guarantees the stability of the lightest particle. Furthermore this 
symmetry permits only annihilation of pairs of DM particles into SM particles. Second, 



the annihilation cross section of DM particles is taken to be typically of electroweak 
strength such that DM particles are in thermal equilibrium in the early Universe and 
eventually freeze-out as the Universe cools down. Third, the contribution of the relativistic 
degrees of freedom to the entropy density and energy density of the Universe are taken 
to go as and respectively. 

In the following we will modify some of these assumptions by introducing the possibility 
of asymmetric DM and of new annihilation processes called semi-annihilation where two 
DM particles annihilate into another DM particle and a SM particle. We will furthermore 
include the possibility of DM annihilation into 3-body processes. We will continue to 
assume that DM annihilation is typically of electroweak strength and keep the case of very 
weak interactions for further development. We will also postpone the full implementation 
of a modification of the entropy and energy contributions to the Universe for a further 
version but will discuss the possibility to modify the temperature dependence of the 
effective number of degrees of freedom. 



2.1 Asymmetric dark matter 

We consider a dark matter particle that is not self-conjugate. This occur for example when 
the DM particle is a Dirac fermion [34] or a complex scalar [35, 36]. The Boltzmann 
equations for the abundance of a DM particle -defined as the number density divided by 
the entropy density - have to be generalized [37, 38, 39]. Recall that for a self-conjugate 
DM particle, X, 

dY^^ <av> , 

ds 3H ^ ' 

where < cif > is the relativistic thermally averaged annihilation cross-section, Y^q is the 
equilibrium abundance, s the entropy and H the Hubble constant. 

We define and Y~ as the abundances of the DM particle/anti-particle {x/x)- The 
evolution equation is then 

with Y^{T) = Y^g{T) = l/2Yeq{T). Note that the factor of 2 in the numerator of 
Eq. 2 appears because we assume that annihilation of pairs of DM(anti-DM) particles 
are zero, thus the thermally average cross section depends only on XX annihilations, 

< av >= \ < av^' >= \ < air+ >. 



Subtracting the two equations in Eq. 2 we obtain 

dY+ dY~ 



(3) 



ds ds 

implying that the difference of abundances AF = Y^ — Y~ is constant. 

To proceed, we define the quantity Y = '1\JY^Y~ which is equivalent to the abundance 
Y = Y^ for a self-conjugate DM. Using the abundance for each species. 



Y± ^ ±Ay + yy^ + 



4 



the evolution equation for Y is derived 



dY < (TV >Y 
ds 




(5) 

This equation is similar to the equation for a self-conjugate DM candidate save for an 
additional term. It can be solved numerically using the usual procedure [40], to ob- 
tain the abundances today, Yq = Y(T = To). The total relic density of the DM parti- 
cle/antiparticle is then given by 



87^ x/Fr? + AY^ 



Qh'' = ^ ° ' (6) 

3-^100 ^Planck So 




where sq is the entropy today, Hiqq = lOOkm/s/Mpc (the Hubble constant today is 
Hq = HiQoh), and Mpianck is the Planck mass. The relic density of each specie can be 
expressed in terms of the asymmetry parameter 6dm, 

n^h' = '"'^ (7) 

with 

6dm - 

Note that the asymmetry, AF, will always increase the relic abundance as compared to 
the one obtained in the same model but without an asymmetry. 

To illustrate the impact of the asymmetry on the relic density, we use a simple ex- 
tension of the SM with one extra lepton generation, 64, 2/4 where z/4 is a massive Dirac 
neutrino. These new leptons interact with SM gauge bosons but do not mix with the 
first three generations. We assume an ad-hoc discrete symmetry where only 64, 2/4 have 
i? = — 1 while SM particles have R = \. Thus when z/4 is lighter than 64 it is stable 
and a potential DM candidate. The dependence of on AF for several DM masses is 
displayed in Fig. 1. ^ When AF is small there is no changes in the symmetric component 
and the asymmetry will lead to a small increase in the relic abundance. When the asym- 
metry times av is large, F_ is very small and the relic density is determined by AF, this 
shows up as a linear dependence of the relic density with AF in Fig. 1. Note that since in 
this case one DM component has disappeared completely, a signal in indirect detection is 
not expected. The linear dependence starts at lower values of AF for larger DM masses, 
since in this model av also increases with the DM mass up to vhdm ~ 50 GeV. Note 
that when finding the starting point for the solution of the evolution equation we assume 
that F — Feg ^ Fgg, when AF times av is very large this assumption will fail, this means 
that the freeze-out picture is inadequate. In this example av is strongly enhanced when 
"^DM ~ ^z/'^ and our solution method breaks when AF 10~^° for this mass. 



^This model is only illustrative and would need to include new quarks to be anomaly free, furthermore 
the light DM masses are incompatible with the bounds on the Higgs invisible width that can be obtained 
from LHC results [41, 42]. 




Figure 1: Qh"^ as a function of AY for different DM masses in the 4th generation Dirac 
neutrino dark matter model for DM masses of 3, 6, 10, 50, 100 GeV as indicated. 



In micrOMEGAs the asymmetry AY is a global parameter, deltaY, that can be set 
anywhere in the code. The relic density computation takes it into account automatically 
to obtain both Qh'^ and Sdm, Eq- 8. The latter is also a global parameter DMasymm, 
see Table (2). Furthermore the parameter DMasymm will have a direct impact on other 
observables computed in micrOMEGAs. 

The computation of the rate for dark matter scattering on nuclei is normally done 
averaging the cross section for DM and anti-DM scattering on nucleons and depends on 
the local DM density which is the sum of DM [p^) and anti-DM (p^)density. When 
the asymmetry parameter Sum is non-zero, the rate will take into account the relative 
abundance of matter and anti-matter given the value of the local DM density specified 
as a global parameter. For indirect detection, the production rate of particles from DM 
annihilation is also modified to take into account the density of DM and anti-DM, 

Q = Uav}^'^, (9) 
2 ah 

and the local density is + p^. Note that the parameter DMasymm=5/)A/, which charac- 
terizes the asymmetry between the relic density of the two components, is also a global 
parameter. If reset after the relic density computation, it will be taken into account for 
computing direct and indirect detection rates and thus will allow the user to examine the 
impact of an asymmetry on these observables. 



2.2 Semi-annihilation 

With discrete symmetries larger than Z2, one introduces new DM annihilation processes as 
well as the possibility of having more than one DM candidate corresponding to the lightest 
particle of two distinct dark sectors. Although micrOMEGAs has been modified to include 
more than one DM candidate [28], this possibility has been implemented only within the 
context of some specific DM models and is not included in this public distribution. In 
this section we therefore discuss only the impact of new processes of the type XX ~^ 



where x is any particle from the dark sector (x correspond to anti-particles) and X is any 
SM particle. These processes are called semi-annihilation [43, 44, 28] 

The simplest example of a discrete symmetry which leads to semi-annihilation pro- 
cesses is Z^. The SM particles have Z3 charge 0, while particles (anti-particles) in the 
dark sector have Z3 charges 1 and —1 respectively. Hence such models feature as usual 
only one DM candidate. With semi-annihilations, the equation for the number density 
reads 

^ = _/^^XX^^X^ / 2 _ 2X _ i/^^XX^xX\ ^^2 _ ) _ ^ffn. (10) 

at ^ ^ ^ 2 ^ 

We define 

a, = + -(w^"^*^) and a = (11) 

2 2 (Tt, 

which means that < a < 1. In terms of the abundance, Y = n/s, where s is the entropy 
density, we obtain 

«)i;q') (12) 

«) V) • (13) 

To solve this equation we follow the usual procedure [40]. Writing Y = Y^q + 6Y we 
find the starting point for the numerical solution of this equation with the Runge-Kutta 
method using 

dY 

3H^ = a,Y,q5Yi2-a), (14) 

where 6Y <^ Y. This is similar to the standard case except that 6Y -> 6Y{1 — a/2). 
Furthermore, when solving numerically the evolution equation, the decoupling condition 
> Fcq^ is modified to 

Y^:^aYY,q + il-a)Y,q\ (15) 

This implies that the freeze-out starts at an earlier time and lasts until a later time 
as compared with the standard case. Although semi-annihilation processes can play a 
significant role in the computation of the relic density, the solution for the abundance 
depends only weakly on the parameter a, typically only by a few percent. This means 
in particular that the standard freeze-out approximation works with a good precision. 
Explicit examples of the impact of semi-annihilation can be found in [28, 36]. The 
function onechannel (see section 2.5) can be used to obtain the contribution of each 
semi-annihilation channel and from there extract a. 

There is no special need to indicate to micrOMEGAs that a model contains semi- 
annihilation processes, as long as all particles in the dark sector are properly identified, 
micrOMEGAs will find the relevant processes. Note that the modified Boltzmann equation 
including both asymmetric DM and semi-annihilation is implemented in the code, however 
both terms cannot be present simultaneously in a specific model. 



or 



— = -sa, {Y^ - aFFeq - (1 
dY 

3H— = a, {Y^ - aFFeq - (1 ■ 



2.3 Thermodynamics of the Universe 

The formulas that give the contribution of one relativistic degree of freedom to the entropy 
density and energy density of the Universe are rather simple. For bosons. 
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while he corresponding functions for fermions have an additional factor |. In the generic 
case one can write 

P{T) = ^T%jfiT) , siT) = -^T'KffiT) (17) 

where geff{T)/heff(T) count the effective numbers of degrees of freedom.^ At freeze-out 
temperatures, the contribution of DM should be strongly suppressed, only SM particles 
contribute to Qeff and hef/-^ At low temperature (T <C lOOMeV), the effective number 
of degree of freedom is defined by photons, neutrinos, light leptons. As the temperature 
increases, one has to add the quarks degrees of freedom. There is a problem of correct 
matching between the low and high temperature region. This issue was studied in [45] 
where other functions which feature different matching between the two regions have 
been derived. By default micrOMEGAs uses the tabulated functions for geff(T)/heff(T) 
obtained in [46, 47], these tables assume that Tqcd = 150MeV. The temperature depen- 
dence of heff is displayed in Fig. 2. To obtain g^ffiT) we solve the following equation 
using the Tables for hefj starting at Tq = 6MeV, 



-^^(ft./,m-*„(r))-H-^j (18) 

The tabulated functions of [45] are also available. heff(T) for three of those functions, 
B, B2 and B3, are compared to the default function in Fig. 2, numerically the difference 
is small. A function (loadHef f Gef f ) allows to substitute any table, see section 2.5. The 
different functions for /ig// can induce a shift in the value of the relic density, typically of a 
few percent. The larger shifts are expected for light DM masses, especially for those where 
the freeze-out temperature is near the QCD transition region (that is Tfo ~ ttldm/'^^ ~ 
150MeV). As an example, we consider the MSSM and display in Fig. 2 the resulting shift 
in Qh'^ for different choices of heff{T) as a function of the neutralino mass. The neutralino 
mass is varied by changing Mi while keeping M^a = Mi -|-51GeV, for each point the value 
of fi is adjusted such that flh^ ~ 0.1147 [48]. Other MSSM soft parameters are fixed to 
Ms = mCeV, Ma = 1.8TeV, Mu = Mm = 500GeV,MQi = 2TeV,Mui = = 1.5TeV, 
At = Ab = -2.8TeV, Ma = 700GeV, tan/3 = 35. 



2.4 Three-body processes 

The pair annihilation of DM particles into a pair of SM particles may not be allowed for 
some kinematic and/or symmetry reasons. A known example is the suppression of the 
s-wave production of massless fermion in the case of a Majorana DM. This suppression is 
lifted by the emission of photons in association with the fermion pair. Another example 
to which we now turn is the case where annihilation into gauge bosons is kinematically 
suppressed. The cross section for annihilation into off-shell W's can then be as large 
as that of other 2 body annihilation channels. This effect was shown to be particularly 
important in models such as the inert doublet model where the DM candidate is a scalar 
and has a standard SU(2) coupling to W's [49, 50]. 

^geff{T) and heff(T) are related via the condition ^^^-^ = '^^l^ 

^We assume that at low temperatures when DM contributes significantly to the matter density, the 
formation of DM has been completed and Y{T) = Const. 
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Figure 2: Left panel : Values of h^ff as a function of T for different QCD equation of states, 
corresponding to the default micrOMEGAs table with Tqcd = ISOMeV (black/full), and to 
case B (blue/dot-dash), B2(red/dash), B3(dot/green) of Ref . [45]. Right panel : Relative 
difference between Vth"^ computed with the default table for hf,ff and cases B,B2,B3, same 
color code as left panel. 



We include a switch that allows to include three and four-body final states involving 
one/two off-shell gauge boson (W/Z). A similar procedure is also used for the decay of 
a Higgs into 3/4-body final states. For DM annihilation or co- annihilation, we consider 
processes such as 

XX ^ XV (19) 

where V = W, Z. To take into account the process with an off-shell V, we use CalcHEP 
to compute the 3-body process 

XX' ^ XII' (20) 

where //' is one leptonic decay channel of V. To obtain the total annihilation cross section 
into a virtual V, we divide this result by the branching ratio, Br{V — )■ W). In general for 
/, I' we use leptons of the first generation except when X is a first generation lepton (this 
can happen in some co-annihilation processes). In this case we use second generation 
leptons to avoid additional poles in matrix elements from unwanted diagrams that do not 
contain a virtual V. For processes involving a virtual Z we use //' = uu to avoid extra 
diagrams with virtual photons, even though those are expected to be small due to the 
large virtuality of the photon. 

Processes with three-body final states can contain diagrams that do not involve a 
virtual V, our procedure for reconstructing the total 3-body cross section from only one 
specific final state would in this case be spoiled. Ideally we would need to select only 
the diagrams with a virtual V corresponding to the process in Eq. 19; this however can 
induce breaking of SU{2) gauge invariance. To solve this problem we select diagrams 
where the total number of virtual vector bosons (W and Z) is at least 1, thus ensuring a 
gauge invariant set of diagrams. 

If X is also a SM vector boson then the procedure we have just described does not 
work very well since one or both of the V can be off-shell. First, note that far below the 
threshold for VV production, the total cross section is twice the one for VW since each 
vector boson can be off-shell, while very near threshold the 3-body result should be close 
to the 2-body result. We therefore need to compute the cross section for four-body final 
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states. This is however computer-time consuming and would significantly slow down our 
code to a point where large scans of parameter space cannot be realized. We therefore use 
a trick to evaluate the 4-body process from the 3-body result after applying a K-factor. 
To obtain the K factor we assume that the 4-body matrix element is proportional to 



{mim2 f + ^ (s - ml - ml)' 

o 



((m? - )2 - rfmf ) {{ml - m!^Y - Vlm!^) 



(21) 



where m^^m^ are virtual masses and m'^^m'^ are pole masses, ri,r2 are particle widths 
and s = 4m^. Note that this is the form of the matrix element for Higgs decay into 
4-body final state through two virtual gauge bosons. We integrate this matrix element 
with Fi = in order to simulate the 3-body result, then we integrate the same matrix 
element with the correct value for Fi. The ratio of these integrals gives the K factor. We 
have compared the value of the annihilation cross section at freeze-out obtained using the 
K-factor with the one obtained with the complete 4-body kinematics, in the MSSM for 
the benchmark point described below the difference is less than 5%. 

We consider only the contribution from virtual W's and Z's. The contribution from a 
virtual Higgs is expected to be small because the width is very small. Contributions from 
virtual top quarks could also be important for DM masses below the top threshold [49], 
such contributions will be included in a future upgrade. 

In Fig. 3 we compare the DM relic density computed using 3-body final state with 
the one using only 2-body final states in the MSSM. Here the LSP mass is changed by 
varying Mi, other soft parameters are fixed to M2 = 200GeV, M3 = ITeV, fi = 126GeV, 
Ma = ITeV, tan/3 = 10, Mu = Mm = Mq^ = M^i = M^i = l.STeV, At = 2.8TeV. 
Differences up to 17% are found for masses just below the W-pair production. Furthermore 
differences up to 10% are also found for LSP masses below mh/2. Note however that in 
this example the relic density is very small for masses below mh/2 since the Higgs exchange 
in s-channel gives an important contribution to DM annihilation. When annihilation is 
completely dominated by Higgs exchange, we find basically the same result for the two 
options, even though the Higgs width is quite different in the two cases. This is because 
the annihilation cross section does not depend on the value of the Higgs width. Note that 
in micrOMEGAs both the Higgs width and the DM pair annihilation are computed with the 
same approximation (i.e. the two-body option does not include the off-shell boson final 
states in the Higgs width). Therefore the difference between the 2-body approximation 
and the calculation that includes off-shell W/Z's is not as large as found in Ref. [51]. 



2.5 New functions and switches 

As in previous versions, the relic density computation is performed with the function 
darkOmega as described in [52], new features include the 3-body final state, which can 
be taken into account by specifying 
•VVdecay 

Switch to turn on processes with off-shell gauge bosons in final states (VVdecay=l) for 
DM annihilation and Higgs decays. micrOMEGAs calculates the width of each particle 
only once. A second call to the function pWidth (whether an explicit call or within the 



in 



10 
10 ^ 

10 " 

-A 

10 

50 55 60 65 70 75 80 

MoM(GeV) 

Figure 3: Left panel: VLh? as a function of Mdm in the MSSM (full) and relative difference 
between the 3-body and 2-body value (dash). The soft parameters used are specified in 
the text. 

computation of a cross section) will return the same result even if the user has changed 
the VVdecay switch. There are two ways to clean the decay table and re- activate the 
VVdecay switch. First one can use the command 
•clearDecayTable 

After this command micrOMEGAs will recalculate particle widths according to the new 
value of VVdecay. Second the sortOddParticles command also cleans all decay tables, 
note that this command must be used to recompute the particle spectrum after changing 
the model parameters. 

If particle widths were stored in the SLHA file downloaded by micrOMEGAs, then the 
SLHA value will be used by micrOMEGAs. To prevent downloading of decays one can 
download SLHA file using slhaRead(/i/eA^ame,mode=4) , see description in Section 7. 

Two new global parameters are available to set the DM asymmetry: 
•deltaY 

Global parameter that describes the difference between the DM and anti-DM abundances 
•DMasymm 

Global parameter that characterizes the relic density of the DM and of the anti-DM 
components, if set after the relic density computation will be taken into account for 
direct and indirect detection rates. 

The temperature dependence of the effective number of degrees of freedom can be set 
with 

• loadHef f Gef f (char*f name) 

loads the file fname which contains a table of heff{T),geff{T) . A positive return value 
corresponds to the number of lines in the table. A negative return value indicates the line 
which creates a problem (e.g. wrong format), the routine returns zero when the file fname 
cannot be opened. The default file is std_thg.tab and is downloaded automatically if 
loadHeffGeff is not called, five other files are provided in the sources/data directory from 
[45]. The user can substitute his/her own table as well, if so, the file must contain three 
columns containing the numerical values for T, hgff, Qeff, the data file can also contain 

1 1 




comments on lines starting with #. 

Several new functions that give information on the contribution of different channels 
to DM annihilation have been added in this version 

• vSigma(T,Beps,f ast) 

calculates the thermally averaged cross section for DM annihilation times velocity at a 
temperature T [GeV], see Eq.2.6 in [53]. The value for av is expressed in [pb]. The 
parameter Beps defines the criteria for including coannihilation channels. The fast=l/0 
option switches between the fast/ accurate calculation. The global array vSigmaTCh con- 
tains the contribution of different channels to vSigma. vSigmaTCh [i] .weight specifies 
the relative weight of the i*^ channel 

vSigmaTCh [i] .prtclEj] (j=0, 4) define particles names for the i*^ channel. 
The last record in vSigmaTCh array has zero weight and NULL particle names. In the 
Fortran version, the function vSigmaTCh(i , weight .particles) with i > 1, serves the 
same purpose. This function returns if i exceeds the number of annihilation channels 
and 1 otherwise, real* 8 weight gives the relative contribution of each annihilation channel. 
character*10 particles(4) contains the names of particles in the annihilation process. 

• oneChanneKXf ,Beps,pl,p2,p3,p4) 

calculates the relative contribution of the channel pl,p2 — p3,p4 to {Qh"^)^^ where 
pi, . . . p4 are particle names. To make a summation over several channels one can 
substitute "*" instead of the particle name. 

• omegaCh 

is an array that contains the relative contribution and particle names for each annihilation 
channel. In the Fortran version one uses instead the function 

omegaChd , weight , particles) . These arrays are similar to vSigmaTCh described above. 
The array omegaCh is filled after calling either darkOmegaFO or printChannels. 

3 DM in neutrino telescopes 

After being captured, DM particles concentrate in the center of the Sun/Earth and then 
annihilate into Standard Model particles. These SM particles further decay producing 
neutrinos that can be observed at the Earth. Annihilation of DM particles in the galaxy 
is also a source of neutrinos. As for other indirect searches, the shape of the neutrino flux 
depends on the dominant DM annihilation channel into SM particles. Direct annihilation 
of DM into neutrinos gives the hardest spectrum followed by muons and W's. As we 
will see below the signal for neutrinos coming from DM captured in a celestial body is 
essentially determined by the cross section for DM scattering on nuclei and is thus related 
to direct detection searches while galactic neutrinos probe the DM annihilation cross 
section. 

The neutrino signal from DM captured in the Sun or Earth usually dominates over 
the one originating from DM annihilation in the galaxies. The best limit on galactic 
neutrinos have been achieved by IceCUBE [54], however the limit on the DM annihilation 
cross section are still orders of magnitudes above the canonical cross section required 
to achieve the relic density assuming a standard cosmological scenario. Limits on the 
muon flux from neutrinos captured in the Sun were obtained by Super-Kamiokande, [15], 
Baksan [55], Amanda and ICECUBE [17], and are currently more sensitive than direct 
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detection experiments in probing DM nuclei spin dependent cross sections on protons. 
The best limit on the SD cross section on protons is currently obtained by ICECUBE, 
(T^^ = 1.2 X 10-^°cm-2 for moM = 200GeV [17]. Note that ICECUBE with the DeepCore 
extension has a threshold energy of lOGeV [56] while Super-Kamiokande has a threshold 
of 1.6 GeV and is thus the best detector for light DM candidates. The neutrino flux and 
the induced muon flux from DM captured in the Sun and the Earth is a new module and 
will be described below. 



3.1 Capture rate 

The capture rate for DM particles in the core of the Sun/Earth depends on the DM- 
nucleus scattering cross-section, as well as on the DM velocity distribution and local 
density. A DM particle will get captured if after scattering with a nucleus at a dis- 
tance r from the center, its velocity is lower than the escape velocity, Vesd^), defined via 
gravitational potential 

= <P{r) (22) 

The capture rate was derived in [57] 

where u is the WIMP velocity at infinity, j[u) is the WIMP velocity distribution normal- 
ized by the condition 

/■oo 

uf{u)du = 1 (24) 



a^A is the scattering cross section for DM on nuclei, ua is the number density of nucleus 
A, and 

aA--m^mARA, Pa- ^^^^^ (25) 

where rriA is the mass and Ra the radius of nucleus A, The computation of a^A takes 
into account the fact that the WIMP scattering cross section is velocity dependent which 
implies a form factor dependence 



\F{q')\' = expi^-Y^j ^^^^ 

where q is the momentum transfer. For the Gaussian form factor^ we assume that 

Ra = (0.91(m^/GeV)^/^ + 0.3) x lO^^^cm (27) 

For spin-independent interactions which add coherently, the sum is performed over 
all nuclei up to ^^Ni. For spin- dependent interaction, there is no coherence effect oc A^, 
therefore hydrogen gives the largest contribution, in this case we do not include the form 
factor. The capture rate can be modified by finite temperature effects, however this effect 
is below 1% [57] so we neglect it. To calculate UA^r) we use the profiles given [58] and 
[59]. 



^Note that in the direct detection module we use a Woods Saxon form factor rather than a Gaussian 
form factor as it is more accurate especially at large , however the Gaussian form factor can be integrated 
analytically. 



3.2 Annihilation 



Assuming a Boltzmann distribution of DM captured by the Sun/Earth, the DM annihi- 
/"^A^ where N^^ 

{av) 



lation rate is N^Ay where Ny is the number of DM particles and 



^eff 

where Veff is the effective volume of the Sun/Earth defined by equation 



Veff = J r^dre J / yin J r^dre j (29) 

and 0(r) is the gravitational potential. The WIMPs are trapped close to the center of the 
Sun, in this region we can consider that the temperature and the density are constant. 
We choose = T(f) where f is the mean radius of the WIMP distribution 

6T(f) ^ 



with G the Newton's constant. The mean radius can be solved numerically from the 
temperature and density profile of the Sun/Earth. 

If we assume that the matter density is constant near the centre of the Sun, the 
gravitational potential simplifies to (f) = 27rGp(f)r^/3, and the effective volume can be 
integrated symbolically to give roughly [60] 

V®^ = 5.8 X 10^°cm=^(m^/GeV)~^/2 V®^ = 2.1 x lO'^'^cm^mJGeV)-'^/^ (31) 
3.3 Evaporation 

DM captured in the Sun/Earth can scatter with nuclei inside the Sun/Earth and gain 
enough energy to escape from the Sun/Earth. Evaporation depends very sensitively on the 
DM mass, it is important for light DM. This effect was estimated in [61, 60], approximate 
formulas for the evaporation rate were derived in [62]. The evaporation rate [62] 

= ^.!^e-x^(^^)/-(^-) (32) 
Tj-'^r-^ I [r) 

where ^^^(j) is the DM escape energy from the center of the Sun/Earth, f is the mean 
radius (Eq. 30) and v the mean speed of DM, 

V = 33 
V ™x / 

T(r) is the temperature profile of the Sun/Earth as given in the standard solar model [58]. 
The evaporation rate depends also on the integrated scattering cross section of DM on 
nuclei inside a radius Tc, ctc, 

pre 

/ V a^AriA^irr'^dr (34) 
Jo A 

where is taken to be the radius where the temperature drops to 95% of its value, 
T{rc) = 0.95T(f). Typically, evaporation in the Sun affects significantly DM particles 
lighter than 3 GeV and is irrelevant for heavier DM particles. 



(Jr. 
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3.4 Neutrino flux 



When the DM is self-conjugate, the equation describing the evolution of the number of 
DM particles A^^ is simply 

4 = C^-A^xK-^xNx, (35) 

where the capture rate, the annihilation rate and the evaporation rate were given above, 
Eqs. 23,28,32. When the evaporation is negligible and the capture and annihilation rates 
are sufficiently large, equilibrium is reached and the annihilation rate is only determined 
by the capture rate. 

Txx = ^^xx^^x = (36) 



The condition for equilibrium is ^y C^A^^t ^ 1. For the Sun, where ^0 = 4.57 x lO^yrs, 
equilibrium is easily reached, for the Earth, equilibrium might not be reached. To take 
into account the most general case, we do not make assumptions about equilibrium and 
solve Eq. 35 numerically. 

In models where the DM is not self-conjugate, one can get different capture rates for 
particles and antiparticles. Furthermore both particle-particle (antiparticle-antiparticle) 
and particle-antiparticle annihilation channels exist. For example for a Dirac fermion 
DM one can get XX ~^ i^xx) ^^'^ XX ~^ (^xx) "where X is any SM particle. The 
equations describing the evolution of the number of DM (anti-)particles, Ny^[N^), are then 

^x = C^-^xxN^N^-^xxNl-EN^, (37) 

The capture and annihilation rates can be evaluated as above, and the two coupled equa- 
tions solved numerically to extract the number density of DM and anti-DM. In the case 
that evaporation is negligible and equilibrium is reached, these equations can be solved 
analytically, see Appendix B. 

After solving for the number density, one can compute the neutrino flux at the Earth. 
First define, 

■'^xx ~ ' ^xx ~ 2^^^'^x ' ^xx ~ ^xx-^x-^x (38) 

The total neutrino spectrum at the Earth, assuming self-annihilation channels are 
solely into neutrino pairs, is given by 

dE, ~ 4nd' I ^ ^^^"^-^^ dE + ^ « ^"^ff dE 



f 



where d = d^ = 1.5 x 10® km for the Sun and d = d^ = 6.378 km for the Earth, 
Br^y is the branching fraction for annihilation into neutrino pairs, Br^j the branching 
fraction into each particle/antiparticle final state //. A'^^ and Nyy{Npp) are the neutrino 
spectra resulting from those annihilations. The neutrino spectra originating from different 

1.K 



annihilation channels into SM particles and taking into account oscillations and Sun/Earth 
medium effects were computed in [63], we use the tables given there. Note that for the 
neutrino pair production an average over the three flavors in the annihilation process is 
assumed. In models where light DM annihilate directly into neutrinos of a given flavor 
(e.g. the light sneutrino scenario in [64]) this is not the case, however this is still a good 
approximation since almost perfect 3-generation mixing is expected for neutrinos below 
10 GeV propagating in the Sun [63]. 



3.5 Muon fluxes 

Neutrinos that reach the Earth interact with rock below or water/ice within the detector 
and will give muons. To compare with the data, one must therefore compute the muon 
flux. First consider muons that are created inside the detector, leading to the contained 
muon flux, [65] 

PmN, ! dE.^^-^ (40) 



dE^ "dE, dE^ 

where pm is the water density, the Avogadro number, and o"cc the cross section for 
charged current interactions for neutrinos on nucleons, vN — )■ pN' . This cross section is 
averaged over protons and neutrons. 

When neutrinos interact with matter outside the detector, one must also take into 
account the fact that muons lose energy on their way to the detector. The energy loss for 
a muon is assumed to be 

dE 

^ = -p,,(a + ^E,) (41) 

where pm = 2.6 g/cm^ is the rock density and a = 2x 10^^ GeVcm^/g, /3 = 3x 10~^ cm^/g 
characterize the average energy loss of the muon traveling through rock or water. A muon 
of energy E'^ will, after propagating a distance z have an energy E"^, where 

and its survival probability will be 



where y = m^/ {crapM) and r is the muon lifetime. 
The upward muon flux is then [65] 



dE^ 



PmNj 



j ^^'^^ / dz'-^^e^'^P..^.{E'^^E,) (44) 



where 

PPm \a + E^f3 



>og I ^ ) (45) 



Note that E'^{z = Zmax) = E^. 

The event rate takes into account the effective area of the detector. For Super-K, the 
detector is cyhndrical with a radius R = 18.9 m and a height 36.2 m. The muons that are 
stopped in the first 7 m are not observed, this corresponds to the muon energy threshold 
of E = 1.6 GeV. Furthermore, only muons of energy larger than 7.7 GeV go through the 
detector. This means that for models with DM masses below 8 GeV, almost all the muons 
will be stopped within the detector. 

Because of oscillations, the neutrino flux at the Earth will contain an almost equal 
mixture of electron, muon and tau neutrinos. Muon neutrinos will convert into muons 
detectable in a detector while electron neutrinos will produce electrons which will not 
propagate for very long. Finally tau neutrinos will produce taus which rapidly decay, 
one of the tau decay modes involves muons and will therefore induce a correction to the 
i/^-induced muon flux. We do not take this effect into account. 

3.6 Functions for neutrino flux 

• neutrinoFlux(f ,f orSun,nu, nu_bar) 

calculates the muon neutrino/anti- neutrino fluxes near the surface of the Earth. This 
function a) calculates capture, annihilation, and evaporation rates, Eqs. 23,28,32; b) 
solves numerically Eq. 35 or Eq. 37; c) calculates all branchings of DM annihilation 
and substitute them in Eq. 39 to get the fluxes. 

Here f (v) is the DM velocity distribution normalized as in Eq. 24. The units are 
km/s for v and s^/km"^ for f(v). At first approximation, one can use the same Maxwell 
function introduced for direct detection. If f orSun==0 then the flux of neutrinos from the 
Earth is calculated, otherwise this function computes the flux of neutrinos from the Sun 
as described in section 3.1. The calculated fluxes are stored in nu and nu_bar arrays of 
dimension NZ=250. The neutrino fluxes are expressed in [l/Year/km^] . The function 
SpectdNdE(E,nu) returns the differential flux of neutrinos in [l/Year/km^/GeV] , and 

displaySpectrum(nu, "nu from Sun [l/Year/km~2/GeV] " ,Emin,Emax, 1) 
allows to display the corresponding spectrum on the screen. 

One can check the neutrino spectra [63] implemented in micrOMEGAs calling the 
function 

• basicNuSpectra(f orSun,pdgN, outN,nu) 

here pdgN is the PDG number of the particle produced initially in DM annihilation. outN 
should be chosen 1 for muon neutrino and -1 for anti-neutrino, nu is neutrino spectrum 
which can be tested by SpectdNdE(E,nu) function. 

• muonUpward(nu,Nu,rho, muon) 

calculates the muon flux which results from interactions of neutrinos with rocks below the 
detector Eq.(44). Here nu and Nu are input arrays containing the neutrino/anti-neutrino 
fluxes calculated by neutrinoFlux. rho is the Earth density ~ 2.6g/cm^. muon is an 
array which stores the resulting sum of /i"*", fi~ fluxes. SpectdNdE(E,muon) gives the 
differential muon flux in [1/Year/km^/GeV] units. 

• muonContained(nu,Nu,rho , muon) calculates the flux of muons produced in a given 
detector volume, Eq.(40). This function has the same parameters as muonUpward except 
that the outgoing array gives the differential muon flux resulting from neutrinos converted 
to muons in a km^ volume given in [1 /Year/km^ /GeV] units, rho is the density of the 
detector in g/cm^. 
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As in other modules of micrOMEGAs we provide tools for model independent calcula- 
tion of neutrino signal from DM captured. The function 

• captureAux(f ,f orSun, csip, csin, csDp, csDn) 

calculates the number of DM particles captured per second assuming the cross sections 
for spin-independent and spin-dependent interactions with protons and neutrons csIp, 
csIn, csDp, csDn are given as input parameters (in [pb]). A negative value for one of 
the cross sections is interpreted as a destructive interference between the proton and neu- 
tron amplitudes. The first two parameters have the same meaning as in the neutrinoFlux 
routine. The result depends implicitly on the global parameters rhoDM and Mcdm in Table 
2. The annihilation rate for this model independent analysis is obtained from the cap- 
ture rate using Eq. 36, and basicNuSpectra described above can be used to obtain the 
neutrino fluxes assuming a value for each channel branching ratio in Eq.(39). 

Two functions allow to estimate the background from atmospheric neutrinos creating 
muons after interaction with rocks below the detector or with water inside the detector. 

• ATMmuonUpward(cosFi ,E) calculates the sum of muon and anti-muon fluxes that result 
from interaction of atmospheric neutrinos with rocks. Units are [l/Year/km^/GeV/Sr] . 
cosFi is the angle between the direction of observation and the center of the Earth. E is 
the muon energy in GeV. 

• ATMmuonContained(cosFi , E, rho) calculates the flux of muons caused by atmo- 
spheric neutrinos which are produced in a given detector volume in [ 1 /Year/km^ /GeV] . 
rho is the density of the detector in [g/cm'^] . cosFi and E are the same as above. 

4 Dark matter detection 

The modules to compute the direct and indirect detection rates are described in Refs. [66, 
67]. Here we just mention the few upgrades including improved values for the quark 
coefficient in nucleons relevant for direct detection, a new module to compute loop-induced 
annihilation of neutralinos into gauge bosons in the NMSSM and CPVMSSM and new 
functions for introducing DM clumps. 

4.1 Gamma- Ray line 

The annihilation of a pair of dark matter particles into photon pairs, or 7Z gives a 
very distinctive DM signal, a monochromatic line at an energy E = m£,]^[ for 77 or 
E = 'mDM{^ — 'm'z/4:m'j^]^.j) for 7Z. Possible evidence for such lines in FermiLAT data have 
been reported [10] although the result is not yet confirmed by the FermiLAT collaboration 
and the DM origin of these lines remains to be established. Such processes occur at the 
one-loop level and are therefore not computed in CalcHEP for generic models. However a 
module for computing loop-induced neutralino pair annihilation processes in the MSSM 
{XX ~^ 77) lZ ) was provided in previous micrOMEGAs versions [68]. The code to compute 
the amplitudes was obtained with SloopS, an automatic code for one-loop calculations 
in the context of the SM and the MSSM [69, 70]. In this version, a similar module is 
provided for the NMSSM [71] and the CPVMSSM using SloopS. The resulting routine, 
lGamma.exe, reads the [N/CPV]MSSM parameters from an SLHA file, and computes the 
cross sections for xX ~^ 77; 7^- 

The calculation of xX 11 y 1^ is based on the same functions for all MSSM-like 
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models, 

• int loopGamma (double *vcs_gz, double *vcs_gg) 

If the call was successful this function returns zero and assigns to the variables vcs_gz 

and vcs_gg the cross sections av for 7Z and 7,7 in [cim? / s\. The function gammaFlux 

can then be used to calculate the photon fluxes, see example in main . c [F] . Note that to 

activate the calculation of loopGamma one has to uncomment the statement 

#define LoopGAMMA 

at the top of the file main . c [F] . 



4.2 Dark matter profile and clumps 

The indirect DM detection signals depend on the DM density in our Galaxy. The DM 
density is given as the product of the local density at the Sun with the halo profile function 

p(r) = p^FhaUr) (46) 

In micrOMEGAs Pq is a global parameter phoDM and Zhao profile [72] 

F..,.(r)=(^)'(^2±^=)'^ (47) 

with parameters a = 1, /3 = 3, 7 = 1, rc = 20[/cpc] is used by default. The parameters of 
Zhao profile can be reset by 

• setProf ileZhao(a,/3,7,rc) 

The functions to set another density profile reads 

• setHaloProf ile(F/iaZo(r)) 

where Fhaio{f^) is any function from distance from Galactic center, r has to be defined in 
[kpc] units. For instance, setHaloProf ile(hProf ileEinasto) sets Einasto profile 



Fhaio{r) = exp 



r 



where by default a = 0.17, but can be changed by 

• setProf ileEinasto(a) 

The command setHaloProf ile(hProfileZhao) sets back Zhao profile. Note that both 
setProf ileZhao and setProf ileEinasto call setHaloProf ile to install the correspond- 
ing profile. 

Dark matter annihilation in the Galaxy depends on the average of the square of the 
DM density, < >. This quantity can be significantly larger than < p >^ when clumps 
of DM are present [73]. In micrOMEGAs , we use a simple model where fd is a constant 
that characterizes the fraction of the total density due to clumps and where all clumps 
occupy the same volume Vd and have a constant density pd- Then assuming clumps do 
not overlap, we get 

< p2 >= p2 ^ /ciPdp. (48) 

This simple description allows to demonstrate the main effect of clumps: far from the 
Galactic center the rate of DM annihilation falls as p(r) rather than as p(r)^. The pa- 
rameters Pd and fd have zero default values. The routine to change these values is 

• setClumpConst (/ci,pd) 
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To be more general, one could assume that pd and fd depend on the distance from galactic 
center. The effect of clumping is then described by the equation 



< p2 > (r) = p{r){p{r) + pi^„p(r)) (49) 

and the function 

• setRhoClumps(p^f/^p) 

allows to implement a more sophisticated clump structure. 



4.3 Direct detection 

The computation of the direct detection rate was presented in [66]. One important 
theoretical uncertainty in the direct detection rate comes from the quark coefficients in 
the nucleon, and in particular the strange quark content. For SI interactions the operator 
{NlniqipqipqlN) IS interpreted as the contribution of the quark q to the nucleon mss, Mn 
with 

(7V|m,^,^,|iV) = /fM^ (50) 

We define the operators 

cr,rAf = mi{p\uu + dd\p) (51) 
(Xs = ms{p\ss\p) (52) 

where m; = {mu+md)/2. The relative importance of strange quark in nucleon is described 
by the ratio 

y ^ 2(p|ggb) ^ m2a^ ^^^^ 
{p\uu + dd\p) rris a^N 

In earlier micrOMEGAs versions, the strange quark content was extracted using the 
operator ctq = mi{p\uu+dd—2ss\p) , this method leads to large uncertainties in the strange 
quark content and in the direct detection rate [74, 75]. Recent lattice QCD calculations 
have allowed a direct determination of as (see [76]), thus reducing significantly the large 
uncertainty in the estimation of y, furthermore the central value for as is also rather small. 
A compilation of recent lattice results is presented in Table 1. Taking the weighted mean 
of these values {^Xi/af /Y^ '^/^f) we get 

a, = 42 ± 5 MeV and a^^ = 34 ± 2 MeV. (54) 

Considering the better and more direct determination of as it is more convenient to set the 
quark from factor from ag rather than from ctq, other parameters are mu/nid, nis/md, ct^at. 
The quark coefficients are expressed as 

fd = -, f^ = —c^fd Fs=— (55) 



rrid 



1 + ^ ) mp{l + a) "^d. rrip 



where 



^_ {N\uu\N) _ 2z-{z-l)y _ {N\uu\N) - {N\-ss\N) ^ 

{N\dd\N) 2 + {z-l)y {N\dd\N) - {N\ss\N) ' ^ ' 
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A new function to compute the quark coefficients is defined, 
• csLLcSca.la.xQuar'kFF {m^ /rnd,rns /rndjCTT^N , as) 

computes the scalar coefficients for the quark content in the nucleon from the quark mass 
ratios mu/mfi^ms/md as well as from ct^tv and a^. The default values given in Table 3 
in the Appendix are obtained from the central values as = 42MeV, aj^N = 34MeV in 
Eq. 54 and from rriu/md = 0.56, m^/md = 20.2 [77]. Note that the mass ratios estimated 
from lattice results rather than from chiral perturbation theory are somewhat different 
rriu/md = 0.46 ± 0.05, m,/md = 27.5 ± 0.3 [77]. 



Table 1: Nucleon coefficients 



cTs (MeV) 


a^^(MeV) 


Reference 


8 ±21 




[78] 


43± 10 




[79] 


54 ±8 




[80] 


22 ±37 




[78] 


31± 16 


50± 10 


[81] 


34+^^ 


42+f 


[82] 


70 ±68 


31 ±5 


[83] 


22 ±20 


32 ±2 


[84] 


2ltT 


45 ±6 


[85] 


125 ± 59 


43 ±6 


[86] 


50± 18 




[87] 




37± 10 


[88] 



5 The Higgs sector at colliders 

The LEP, Tevatron and LHC colliders have performed extensive searches for the Higgs 
boson(s). Recently the LHC has announced the discovery of a Higgs-like particle with a 
mass m/i = 125.5GeV. [18, 19] A signal was found notably in the two-photon channel and 
in the two gauge bosons decay modes. The observation of this new particle as well as the 
non-observation of a Higgs particle in many other channels at the LHC, LEP or Tevatron 
both put strong constraints on the parameter space of the extensions of the standard 
model. Two new features of micrOMEGAs allow to compute more precisely the partial 
widths of the Higgs: the partial width for the Higgs into off-shell W/Z's are included 
following the approach discussed in section 2.4 and loop- induced couplings of the Higgs 
are included via effective vertices. 

5.1 Loop- induced Higgs production and decay 

The loop-induced couplings of the Higgs to gluons and photons play and important role 
in Higgs production and decay at colliders. Thus the micrOMEGAs models which are 
based on tree-level vertices need to be extended to include these interactions. For this 
we introduce effective operators of the form XhF^j^F^^ for the CP even part and X'hF^y 
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F^'^ for the CP-odd part^. These operators receive contributions from charged scalars, 
vector bosons or fermions. For a generic Lagrangian describing the Higgs interactions 
with fermions ip, scalars and vector bosons V^, 



the contribution to the effective operators read 



(57) 
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(5^ 



(59) 



where q is the electric charge of the virtual particle for 77 operators and is 1 for the gluon 
operator, is the color factor associated with the virtual particle. For the fundamental 
SU(3) representation = 3 for 77 and = 1/2 for gg. For the adjoint representation 
= 8 for photons and = —2 for gluons. The functions A1/2, Ai, Aq, A1/2 are defined 
in [90, 91]. 

There are important QCD corrections to the effective vertices with coloured particles 
in the loop. For h'-ff the QCD corrections are described as an overall factor 



1 



a,{Mh/2) 
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(60) 



where Mi is the mass of the particle in the loop and the Ci functions are known for 
fermions and scalars in the fundamental representation of SU(3). The code HDECAY [92] 
is used to generate tables for these functions. 

The QCD corrections to hgg to order have been computed in the limit of a heavy 
top quark. These include both vertex corrections and corrections from gluon emission [93]. 
The corrections are therefore split into a radiation factor and a vertex factor. For lighter 
quarks or for new colored particles, only NLO vertex corrections are known. To take 
into account QCD corrections for all colored fermions and scalars in the loop, we use the 
following prescription. We use a common radiation factor, R, for all colored particles and 
introduce additional correction factors in the amplitude for top/heavy new quarks (Ct), 
light quarks {Cq) or colored scalars (Cg). The hgg vertex including QCD corrections for 
SM fermions and colored scalars (such as squarks in the MSSM) thus reads 



A 



QCD 



-R 



(61) 



where A\^yr denote the contribution of particle X to Eq. 58. The QCD factors with 
a = as(Q)/vr, read [93], 



R^ = 1 



149 
~12 



a + 68.6480^ - 212.447a^ 



(62) 



^These operators arc also used in CalCHEP 
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and 



= 1 + — a + (6. 1537 - 2.8542 log (^^j)a^ 
+ (10.999 - 17.93 log -| + 5.47 log 




(63) 



= 1 + — a 
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(64) 




= 1 + -a, 
2 ' 



(65) 



where Q = Mh/2. The top contribution to the hgg partial widths agrees with Ref. [93] 
to order a^. Our approach which consists in a modification of the effective /177 and 
hgg vertices, necessarily incorporates only part of the higher order corrections to the 
partial widths (for example for lighter quarks the contribution to order are taken into 
account while only some of the contributions are included). Furthermore with the 
effective vertex approach we cannot introduce QCD corrections for the interference term 
between two contributions. It is therefore expected that we will have differences with 
other treatment of QCD corrections. We indeed observe some differences with HDECAY. 
For instance for lighter quarks or for new particles, the QCD factor in the partial width 
is in agreement with HDECAY only to order a^. Nevertheless for a SM Higgs of 126 GeV, 
differences with HDECAY are below the percent level. For the SM-like Higgs of the MSSM, 
differences are also generally small (at the few percent level). Larger discrepancies can 
be found for the heavier Higgs of the MSSM, in particular when the parameters are 
chosen such that there is a strong interference between the top and stop contribution. In 
supersymmetric models, higher-order SUSY corrections, the so-called Arrib corrections, 
are taken into account by a redefinition of the hbb vertex. These are also included in 
micrOMEGAs when defining the b-quark contribution to h'-f-f and hgg, such corrections 
are not included in the loop-induced processes in HDECAY thus leading to discrepancies 
between the two codes when the b-quark contribution is important. 

In most extensions of the standard model, the Higgs sector contains several doublets 
and/or singlets. For all Higgs scalars and pseudoscalars, we follow the same procedure 
for computing the H'-fy and Hgg vertices, including QCD corrections. Note however 
that these corrections were obtained in the heavy top quark limit which would not be 
appropriate for the Higgses that are very heavy. Nevertheless we use these results for all 
masses, since a high precision is not essential for the loop-induced processes when the 
Higgs is very heavy, indeed both the gluon fusion Higgs production and the decay in the 
77 mode are suppressed at large masses. 

The functions that can be used to compute the effective Higgs vertices are 



to compute the fermion/scalar /vector loop contribution to the hgg,h'j'j vertices, Q = 
{Mh/2Mf where M is the mass of the fermion/scalar /vector in the loop. 



to compute the fermion loop contribution to the Agg, A77 vertices, Q is defined as above. 
• HgamlF(Q) ,HgamlS(Q) 

QCD corrections to the fermion/scalar contribution to the h'j'j effective vertex, see Eq. 60. 



•HggF(Q),HggS(q),HggV(q) 



• HggA(Q) 



• HgamlACQ) 

QCD corrections to the fermion contribution to the ^477 effective vertex, see Eq. 60. 

Implementation of the loop-induced vertices can be realised using LanHEP, the nec- 
essary functions are described next. 

5.2 Implementation of effective vertices 

A new feature in LanHEP allows to extract the expression from any generated vertex. 
This can be used to implement the different contributions to the loop-induced H'-j'-j and 
Hgg vertices. The function which allows to extract a specific coupling is 
Coef Vrt ( [ particles ] , [ options ] ) , 

where particles is the comma-separated list of particles in the vertex, and options specifies 
the tensor structure of the required coefficient as well as other options. The only tensor 
structure that can be extracted are 1, 7^, 7^^, 7^7^,^'^ and the function Coef Vrt must be 
called separately for each tensor structure. The second list can contain the following 
symbols in any order: 

• gamma corresponds to only 7^^ matrices in the vertex, gammaS to only 7^ matrices 
and gamma, gammaS to product of 7^7^ matrices. If none of these symbols are listed, 
then only terms which do not include 7-matrices will be extracted; 

• moment ( P ) where P is a particle name selects the terms depending on this par- 
ticle momentum, note that vertices corresponding to higher order operators and 
containing the product of several momenta cannot be handled; 

• re, im - take real or imaginary part of the expression. It can be combined with 
options above. 

• abbr - use abbreviations in the returned expressions, this option is useful for long 
expressions. 

The following example illustrates how LanHEP can be used to compute the sleptons 
contribution to the /?,—)■ 77 effective vertex in supersymmetry. First, we define terms 
depending on the loop functions : 

_p=[eL,eR,mL,mR,ll,12] in parameter rhS__p=creal(HggS((Q/2/MS_p)**2)) , 
ihS__p=cimag (HggS ( (Q/2/MS_p) **2) ) . 

Here eL,...12 are the names for the three generations of charged sleptons, MSeL...MS12 
are their masses, and we declare the new parameters rhS_eL...ihS_12. The function HggS 
describes the loop integrals in Eq. 58. Then we obtain expressions from the vertices of the 
Higgs boson interaction with these sleptons (divided by the square of the slepton mass): 

_p=[eL,eR,mL,mR,ll,12] in parameter AhS__p=CoefVrt( [anti(_p) ,_p,h] ) 
/(mass _p)**2/2. 

Note that in general it is preferable to specify whether one wants to extract the real or 
imaginary part of the vertex, this option was omitted here because we know that this 
vertex is real. Similarly the option abbr is useful to get a compact expression for the 
vertices involving staus. The charge of particles can, for example, be defined as 
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_p in [eL,eR,mL,mR,ll,12] alias Q(_p)=-1 

Finally, the contribution to the Higgs decay amplitude reads: 

LAAhs=cabs ( (rhS_#sl+i*ihS_#sl) *AhS_#sl*Q (si) *Q (si) 
where sl= [eL, eR,niL,mR, 11 , 12] ) . 

In this statement the inner where keyword is used to make the sum of expressions where 
the list of symbols on the right-hand side is substituted into the template expression on 
the left-hand side. The # sign means concatenation of the symbols to the left-hand and 
right-hand side. To this expression, one must in general add the relevant color factor that 
appear in Eq. 58 as well as the QCD factors in Eq. 61. 



5.3 Higgs couplings and Higgs exclusion 

In addition to a measurement of the Higgs mass, the LHC collaborations provide estimates 
of the signal strength (/i) as compared to the SM expectation in various channels. For 
example for the two-photon decay of a Higgs produced via gluon fusion(ggF) or vector 
boson fusion(VBF) 

ggF ^ (^{99 h)NpBr{h 'y'j)NP vbf ^ (y{WW h)NpBr{h 'J'j)np 

where NP stands for the new physics model under consideration. The ratio of the Higgs 
partial widths can be obtained easily from the coupling squared normalized to the SM 
coupling 

Chxx = glxANP)/9lxx{SM) (67) 
and similarly for loop-induced processes, 

CHGG = \Igg{NP)/\Igg{SM) (68) 

where G = g,'-f and XhCG are given in Eqs. 58,60,61. We make the approximation that the 
ratio of the production cross section is given by the ratio of partial widths even though 
the QCD correction factors can differ in both cases. For example, the signal strength is 
thus given by 

(69) 

J- NP 

The reduced couplings as well as all Higgs partial widths are computed within micrOMEGAs and 
are provided in a SLHA file, HB . slha. The total width for the SM, necessary to compute 
Higgs branching ratios, can be computed from reduced couplings and the total width in 
the NP model. 

To obtain the limits on the Higgs sector for models with more than one Higgs bo- 
son, the results for Higgs signal strengths in different channels described above can be 
compared to the latest results of the LHC. Alternatively an interface to the public code 
HiggsBounds [94] which allows to take into account the exclusion limits provided by the 
experimental collaborations for many different channels is provided. In micrOMEGAs, the 
interface to HiggsBounds is realized via SLHA files, using the effective coupling option for 
specifying the input. More specifically the file HB . slha contains the normalized couplings 
squared of the Higgs to all SM particles, including the normalized couplings squared to 
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77, 7-^) 99 ^ ^ well all Higgs partial widths and the top decay width. A reference value 
is defined for the couplings that do not exist in the standard model, for example the 
reference value for the Z couplings to two Higgses, 9hihz — (c/2s^c^)^. In the MSSM, 
analytical formulas for the effective couplings normalized to the SM ones are used for 
h,H,A — > f f,WW, ZZ while loop-induced couplings, h,H,A — )■ ■y'y,gg, are computed 
using the formulas given above. The same procedure is used for most of the models 
we provide. In the CPVMSSM , the only difference is that loop-induced couplings are 
computed with CPSuperH. In the NMSSM, we have two options for the Higgs sector : 
collider constraints can be computed with NMSSMTools or with HiggsBounds. As for 
other models, the normalized Higgs couplings squared are provided in the Hb.slha file. 

The SLHA file is then fed to HiggsBounds and the output of HiggsBounds is read 
using the slhaVal function [33]. Note that the output of HiggsBounds consist of two 
numbers, HBresult and obsratio which indicate whether a point has been excluded at 
the 95%C.L. by one of the experimental results considered. 
HBresult obsratio 

> 1 parameter point is excluded 

1 < 1 parameter point is not excluded 
-1 < invalid parameter point 

Information on the most constraining channel is also given as output. 

In a generic model implemented by the user, the effective Lagrangian for ^177 and 
hgg have to be included in the model file as described in section 5.1. Furthermore if 
one wants to use the HiggsBounds interface, an SLHA file containing the normalized 
couplings squared must also be provided by the user. In models available in micrOMEGAs, 
the following function can be used for this purpose 
• HBblocksC "filename") 

writes at the end of an SLHA file, the new Blocks that contain the effective couplings of 
the Higgses normalized to the SM ones. 

6 Dark matter models 

Several models are included in the micrOMEGAs distribution: the MSSM and some of its 
extensions (NMSSM and CPVMSSM) as well as the little Higgs model (LHM), a model 
with a right-handed neutrino (RHNM), the inert doublet model (IDM), a model with an 
extended scalar sector and a Z3 symmetry(Z3ID) [28] as well as a model with a fourth 
generation of leptons (SM4). All but the last three models were available in previous 
versions. All models have been extended to include effective higg and /ii77 vertices for all 
Higgses in the model and to output a file containing reduced Higgs couplings (HB.slha) - 
the only exception is the RHNM which features SM Higgs couplings. Both the MSSM and 
NMSSM have been updated, we briefiy describe the main new features of these models. 
A more complete description of all models will appear as a separate publication [95]. 

6.1 MSSM 

The model files were regenerated using SLHAplus [96] and feature an improved Higgs 
sector which includes 

^The coupling Hi jZ will be included in a future release 
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• Effective vertices for Hgg and if 77 where H = h, H, A. 

• Tlireshold corrections for r's in addition to tfiose for b-quarks [97, 98]. 

• Improved Higgs effective potential: as in previous versions we write an effective 
Lagrangian to introduce in a gauge invariant manner the higher-order corrections 
to the Higgs masses and couphngs. The effective scalar potential contains seven 
parameters (Aj) [99], analytic formulas are used to compute Ai, A3, Ag [100, 92] other 
A's are computed from these and the four Higgs masses, mh,mH,TnA,TnH+ given 
by the spectrum calculator. This improvement has an impact on the Higgs self- 
couplings, the effect on the DM relic density is however negligible. 

• The mass limit function checks that the Higgs mass is in the range 123 < nih < 
128 GeV consistent with the LHC results and allowing for a theoretical uncertainty. 
A warning is issued if this is not the case. The Higgs couplings are not explicitly 
checked. 

Other new features of the MSSM model 

• The default spectrum generator has been updated to Suspect_2.4.1 

• In addition to the CMSSM (called SUGRA), AMSB, and EWSB, the user can 
choose to work with non-universal Higgs masses at the GUT scale (SUGRANUH). 
The scalar Higgs masses are not used as input parameters, rather they are replaced 
by /i, defined at the EW symmetry breaking scale and the pseudoscalar pole mass, 
rriA- The input parameters are therefore mo, mi/2, ^o? tan/3, /i, in addition to 
the optional standard model parameters ■mt,mb{mi)),as{'mz)- Defining input pa- 
rameters at different scales follows the standard in SLHA2 [101]. 

• We have updated the routine for 6 — )■ 57 in order to reproduce the SM value at 
NNLO, default values for all parameters are given in the online manual. 

• We have added three flavor observables in the Kaon and D-meson systems, the 
leptonic decays, — )■ fiu and Dg fiiy, ru as described below. Note that these 
constraints are also available via the Superlso interface [102]. 

The two-body kaon leptonic decays receive a contribution from the charged Higgs at 
tree-level. However uncertainties in the Kaon parameters are very large, we therefore 
compute only a quantity that has much smaller uncertainties. This is defined as [103] 



123 



^ m\ f md\ tan/32 



mjj± \ nis J 1 + eo tan (3 



(70) 



with rriK = 0.494 GeV and rrid/ms = 1/20.2. Here the dependence on the supersymmetric 
parameters arises at higher order in perturbation theory (through the eg term) but can 
be large because of the tan/3 enhancement, 

2(1' II m? "^j 
e. = -^H,C\^^) (71) 

y y 

The measured value for this quantity is R123 = 1.004 ± 0.007 [103]. 
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The branching ratios Ds — > tu and Dg — ?■ /iz^ also receive a contribution from the 
charged Higgs, [104] 



r^2 

iTA 12^2 



57r 



nis tan^ (3 



1 



rric 1+eo tan /3 



(72) 



with eo as in Eq. 71 with mj^ — )■ m^^ These branching ratios were measured recently, 
[105] 



Br{Ds — ^iv) 



(5.44 ±0.22) X 10~2 
(5.54 ±0.24) X 10"^ 



(73) 



6.2 NMSSM 

The NMSSM model files have also been updated and a new interface to NMSSMTools_4.0 [29] 
is provided. Thus all new features of NMSSMTools [29, 30] are included, in particular the 
computation of reduced Higgs couplings. The main new features of micrOMEGAs are: 

• Effective vertices for Hgg and H'^'^ where H = hi, aj, i = 1,3; j = 1, 2. 

• Improved Higgs effective potential: as in previous versions we write an effective La- 
grangian to introduce in a gauge invariant manner the higher-order corrections to 
the Higgs masses and couplings. The effective scalar potential contains ten param- 
eters (Aj) [106]. Analytic formulas are used to compute the Aj's which are then fed 
to micrOMEGAs. The new Higgs self-coupling can in some cases have a large effect 
on the DM relic density. 



7 External codes in micrOMEGAs 

The subdirectory Packages in inicrOMEGAs_3.1 contains external public codes such as 
SuSpect [31] , CPsuperH [32] , HiggsBounds-3 .8.1 [94] , LoopTools-2 . 1 [107] . These 
external codes are not compiled during the general installation of micrOMEGAs, rather they 
are compiled either when installing a model requiring one of these codes or in runtime of 
micrOMEGAs executables. The micrOMEGAs interface to these codes (except for Loop- 
Tools) is done via files written following the SLHA convention. Routines of the SLHAplus 
package [96] included in micrOMEGAs can be used to read the information contained in 
these files, as in previous versions the following function is available 

• slhaRead(f ilename ,mode) 

reads all or part of the data from the file filename, mode is an integer which determines 
which part of the data should be read from the file, mode= I*ml+2*m2+4*m4+8*m8+16*ml6 
where 

ml = 0/1 - overwrites all/keeps old data 

m2 = 0/1 - ignore errors in input file/ stop in case of error 

m4 = 0/1 - read DECAY /do not read DECAY 

m8 = 0/1 - read BLOCK/do not read BLOCK 

ml 6 = 0/1 - read QNUMBERS/do not read QNUMBERS 
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For example mode=20 (m4=l ,ml6=l) is an instruction to overwrite all previous data 
and read only the information stored in the BLOCK sections of filename. In the same 
manner mode=25=l+8+16 is an instruction to add information from DECAY to the data 
obtained previously. 

The function slhaRead returns the values: 



- successful reading 

-1 - can not open file 

-2 - invalid data as indicated by SPINFO 

-3 - no data 



n>0 - wrong file format at line n 

In SLHA, information is stored in blocks identifiable with the convention: a line 
starting with the keyword "BLOCK" followed by the name of the block. Each block then 
contains a list of records, each on a separate line started with space symbol, that give 
information on physical parameters. The structure of a block record is not defined in 
SLHA. In the first SLHA implementations, the structure of a block was the following : 
key parameters defining the variables followed by some numerical parameters giving the 
value of these variables. For example, 

BLOCK MASS 

6 173.2 # top quark mass 

1000022 1.34362460e+02 # ~neutralino(l) 

New codes have actually introduced new formats for blocks, for example in the Flavour 
Les Houches Accord [108] key parameters are specified both before and after the numerical 
parameters. 

Block FOBS # Flavour observables 

# ParentPDG type value q NDA IDl 1D2 1D3 ... comment 

5 1 2.83249483e-04 2 3 22 # BR(b->s gamma) 

The SLHAplus package [96] was updated to read blocks with any format, in addition 
to the function slhaVal, we have introduced the function 

• slhaValFormat (char * Block, double Q, char * format) 

which allows to extract numerical information corresponding to the key parameter de- 
scribed in format. The format parameter is actually a C-style format. For instance in 
the "BLOCK MASS" above, the mass of the neutralino can be obtained by 

MNEl=slhaValFormat("MASS", 0., " 1000022 "/.If ") ; 
while in the "BLOCK FOBS", the branching 6 — )■ 57 can be obtained by 
B2sg=slhaValFormat("F0BS" , 0., "5 1 "/.If l*f 2 3 22"); 

The parameter Q is meaningful only for blocks which feature a scale dependence. Note 
that the number of space symbols in format is irrelevant since it has the same properties 
as the C scan/ function. 

A block can also contain a textual information. For example, in HIGGSBOUNDS a block 
contains the following records. 

Block HiggsBoundsResults 

#CHANNELTYPE 1: channel with the highest statistical sensitivity 



1 1 328 # channel id number 

12 1 # HBresult 

1 3 0.72692779334500290 # obsratio 

14 1 # ncombined 

1 5 I I (p p)->h+..., h=l where h is SM-like (CMS-PAS-HIG-12-008) I I # text description of channel 
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In particular, the last record contains the name of the channel which gives the strongest 
constraint on the Higgs. To extract the name of this channel one can use the new function 

• slhaSTRFormatC'HiggsBoundsResults" , "1 5 II % [~ I ] I I " , channel) ; 
which will write the channel name in text parameter channel. 

The following packages have been attached to micrOMEGAs using the SLHAplus inter- 
face: SOFTSUSY [109], SPHENO [110], ISA JET [111], Superlso [102], SUSY-HIT [112], 
HDECAY [92], SusyBSG [113]. 

Note that the make clean command launched in the micrOMEGAs root directory 
subsequently calls make clean for all sub-directories in the Packages directory. 

8 Installation and running the code 

8.1 Installation. 

To download micrOMEGAs , go to 

http : / /lapth . cnrs . f r/micromegas 
and unpack the file received, micromegas_3.1 .tgz, with the command 

tar -xvzf micromegas_3.1 .tgz 
This should create the directory micromegas_3.1/. In case of problems and questions 
micromegasQlapth. cnrs . f r 

To compile CalcHEP and micrOMEGAs go to the micrOMEGAs directory and launch 
gmake 

If gmake is absent, then make should work like gmake. In principle micrOMEGAs de- 
fines automatically the names of C and Fortran compilers and the flags for compila- 
tion. If you meet a problem, open the file which contain the compiler specifications, 
CalcHEP_src/FlagsForSh, improve it, and launch [g]make again. 

After compilation of micrOMEGAs one has to compile the executable to compute DM 
related observables in a specific model. To do this, go to the model directory, say MSSM, 
and launch 

[g]make main=main.c 
It should generate the executable main. In the same manner 

gmake m3iin=filename.ext 
generates the executable filename based on the source file filename. ext. For ext we sup- 
port 3 options: 'c' , 'F', 'cpp' which correspond to C, FORTRAN and C++ sources. [g]make 
called in the model directory automatically launches [g]make in subdirectories lib and 
work to compile 

lib/aLib.a - library of auxiliary model functions, e.g. constraints, 
work/work_aux.a - library of model particles, free and dependent parameters. 

Default versions of mam. c/F programs need some arguments which have to be specified 
in command lines. If launched without arguments main explains which parameter are 
needed. As a rule main needs the name of a file containing the numerical values of the 
free parameters of the model. 

To add a new model to micrOMEGAs one has to launch the command 



./newProject ModelName 
which creates the directory modelName which contains all subdirectories and commands 
required. Then one has to prepare the CalcHEP model files and place them in the sub- 
directory work/models/, furthermore auxiliary codes used in the model have to be in- 
cluded in the sub-directory lib/. 



8.2 Module structure of main program 

Each model included in micrOMEGAs is accompanied with sample files for C and Fortran 
programs which call micrOMEGAs routines, the main.c, main.F G\es^ . These files consist 
of several modules enclosed between the instructions 



#ifdef KEYWORD 



#endif 

where KEYWORD is the name of the module, all modules are completely independent. The 
user can comment or uncomment any set of define instructions to suit his/her need. The 
following modules are available for all models: 

MASSES_INFO //Masses and widths of Higgs particles and masses 
of particles in the odd sector. 
CONSTRAINTS //Checks various constraints. 

HIGGSBOUNDS //Calls HiggsBounds to constrain the Higgs sector. 
OMEGA //DM Relic density. 

1NDIRECT_DETECTI0N //Signals of DM annihilation in the galactic halo. 
RESET.FORMFACTORS //Redefinition of nucleon Form Factors. 
CDM_NUCLEON //Amplitudes and cross sections for DM-nucleon scattering. 
CDM_NUCLEUS //Cross section and recoil energy distribution for 
DM-nucleus scattering. 

NEUTRINO //Neutrino flux from DM captured in the Sun/Earth 

and the corresponding muon flux. 

DECAYS //Decay widths and branching ratios. 

CR0SS_SECT10NS//Cross sections for $2\to2$ processes. 
CLEAN //Cleans intermediate files. 

For the three MSSM-like models, a module to calculate neutralino annihilation in 77 
and 7Z is also available, it is driven by the #def ine LoopGAMMA instruction. 

All models also contain the preprocessor command #def ine SHOWPLOTS which gener- 
ates various plots in CalcHEP. By default it is switched off. 

A complete description of all micrOMEGAs functions is available in the online manual, 
see also [33]. This updated manual is also provided with the code, the file maiiual3.pdf 
can be found in the micrOMEGAs_3.X/man directory. 



^The C++ version of the sample main program can be obtained by copying main.c to main.cpp 



8.3 Sample output 



Running the main.c sample file in the micrOMEGAs/MSSM directory choosing the options 
MASSES.INFO , CONSTRAINTS , OMEGA, INDIRECT.DETECTION, CDM.NUCLEON , NEUTRINO,DECAYS 
and with the mssmh.par input file will lead to the following output. 

========= EWSB scale input ========= 

Spectrum calculator is suspect 
Initial file "mssmh.par" 

The following parapeters keep default values: 

alfSMZ=1.1720E-01 MW=8 . 0200E+01 MZ=9 . 1300E+01 Ml=l . 7770E+00 
McMc=1.2700E+00 MbMb=4 . 2300E+00 Mtp=l . 7350E+02 Au=0 . OOOOE+00 
Ad=0 . OOOOE+00 wt=1.4420E+00 wZ=2 . 5020E+00 wW=2 . 0940E+00 
SU_read_leshouches : end of file 
OUTPUT in SLHA format in suspect2_lha. out 
RUN TERMINATED : OUTPUT in suspect2.out 
Warnings from spectrum calculator: 
none 

Dark matter candidate is '~ol' with spin=l/2 mass=l . 89E+02 

~ol = 0.950*bino -0.062*wino +0 . 257*higgsinol -0 . 167*higgsino2 

=== MASSES OF HIGGS AND SUSY PARTICLES: === 
Higgs masses and widths 

h 126.08 3.49E-03 
H 699.89 5.16E+00 
H3 700.00 5.80E+00 
H+ 704.62 4.74E+00 



Masses of odd sector Particles: 
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==== Physical Constraints : ===== 
deltartho=3 . 56E-04 
gmuon=l . 15E-09 

bsgnlo=3.05E-04 ( SM 3.25E-04 ) 
bsmumu=3 . 38E-09 
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btaunu=9.56E-01 

dtaunu=5 . 17E-02 dmunu=5 . 33E-03 

R123=9.998166E-01 

MassLimits OK 

==== Calculation of relic density ===== 
Xf=2.65e+01 0mega=1.03e-01 

Channels which contribute to 1/ (omega) more than 1%. 
Relative contributions in % are displayed 
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==== Indirect detection ======= 

Channel vcs[cm~3/s] 

annihilation cross section 9.22E-27 cm~3/s 
~ol,~ol -> t T 5.34E-27 
~ol,~ol -> 1 L 1.90E-27 
~ol,~ol -> b B 1.06E-27 
~ol,~ol -> W+ W- 5.57E-28 
~ol,~ol -> Z Z 2.56E-28 
~ol,~ol -> Z h l.lOE-28 
sigmav=9 . 22E-27 [cm"3/s] 

Photon flux for angle of sight f=0.00[rad] 

and spherical region described by cone with angle 0.0349[rad] 
Photon flux = 1.61E-14[cm"2 s GeV] "-[-1} for E=94.5[GeV] 
Positron flux = 5 . 61E-13 [cm"2 sr s GeV] "{-1} for E=94.5[GeV] 
Antiproton flux = 6 . 34E-13 [cm"2 sr s GeV] "{-1} for E=94.5[GeV] 

==== Calculation of CDM-nucleons amplitudes ===== 



CDM-nucleon micrOMEGAs amplitudes: 
proton: SI -3.547E-09 SD -2.115E-07 
neutron: SI -3.634E-09 SD 1.862E-07 
CDM-nucleon cross sections [pb] : 

proton SI 5.446E-09 SD 5.807E-05 

neutron SI 5.715E-09 SD 4.503E-05 

===============Neutrino Telescope======= for Sun 

E>1.0E+00 GeV neutrino flux 5.783577E+10 [l/Year/km~2] 

E>1.0E+00 GeV anti-neutrino flux 6.224851E+10 [l/Year/km"2] 
E>1.0E+00 GeV Upward muon flux 2.258430E+02 [l/Year/km~2] 
E>1.0E+00 GeV Contained muon flux 1.217184E+03 [l/Year/km~3] 

================= Decays ============== 

h : total width=4.564910E-03 
and Branchings : 



2. 


,549952E- 


-02 


h 


-> 


Z,Z 


2. 


,098720E- 


-01 


h 


-> 


w-^,w- 


2. 


. 147891E- 


-03 


h 


-> 


A,A 


6. 


.913403E- 


-02 


h 


-> 


G,G 


6. 


.332254E- 


-02 


h 


-> 


1,L 


1. 


.367769E- 


-04 


h 


-> 


u,U 


1. 


.509459E- 


-04 


h 


-> 


d,D 


2. 


.669341E- 


-02 


h 


-> 


c,C 


1. 


.509459E- 


-04 


h 


-> 


s,S 


6. 


.028919E- 


-01 


h 


-> 


b,B 



~o2 : total width=l . 357505E-01 
and Branchings : 



4, 
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-02 
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z, 


~ol 
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-01 
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~eR 
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1, 
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-01 


~o2 
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m, 


~MR 
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-> 
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1, 


~L1 
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Appendix 

A Analytical solution of neutrino flux 

When DM is not self-conjugate, Eqs. 37 for the number density can be solved analytically 
if we assume that equilibrium is reached and that the evaporation rate is negligible. 
Define 

/3 = ^ ;« = ^ x = ^ (74) 
^xx ^^x 

Then after equilibrium is reached the annihilation rate is determined by the capture rates 
as well as a, here we assume A^^ = A^^, 

1 + ax ^ ^ 

_ C^f xa \ _ X 

XX ~ o \ 1 I 1 ~ X 



2 \1 + ax I X + a 



r - V " 

XX r, X 



where 



X = — 
2a 



2 ""l + ax 



B Global Parameters 

A list of the global parameters and their default values are given in Tables 2, 3. The 
numerical value for any of these parameters can be simply reset anywhere in the code. 



Table 2: Global parameters of micrOMEGAs 



Name 


default value 


units 


comments 


Mcdm 




GeV 


Mass of Dark Matter particle 


deltaY 







Difference of DM/anti-DM abundances 


DMasymm 







Asymmetry between relic density of DM/anti-DM 


K_dif 


0.0112 


kpcVMyr 


The normalized diffusion coefficient 


L_dif 


4 


kpc 


Vertical size of Galaxy diffisive halo 


Delta_dif 


0.7 




Slope of diffision coefficient 


Tau_dif 




s 


Electron energy loss time 


Vc_dif 





km/s 


Convective Galactic vind 


Fermi _a 


0.52 


fm 


nuclei surface thickness 


Fermi_b 


-0.6 


fm 


parameters to set the nuclei radius with 


Fermi _c 


1.23 


fm 


Ra = cA^I^ + h 


Rsun 


8.5 


kpc 


Distance from Sun to center of Galaxy 


Rdisk 


20 


kpc 


Radius of galactic diffusion disk 


rhoDM 


0.3 


GeV/cm^ 


Dark Matter density at Rsun 


Vearth 


225.2 


km/s 


Galaxy velocity of Earth 



Table 3: Global parameters of micrOMEGAs: nucleon quark form factors 



Proton 


Neutron 




Name 


value 


Name 


value 


comments 


ScalarFFPd 


0.0191 


ScalarFFNd 


0.0273 




ScalarFFPu 


0.0153 


ScalarFFNu 


0.011 


Scalar form factor 


ScalarFFPs 


0.0447 


ScalarFFNs 


0.0447 




pVectorFFPd 


-0.427 


pVectorFFNd 


0.842 




pVectorFFPu 


0.842 


pVectorFFNu 


-0.427 


Axial- vector form factor 


pVectorFFPs 


-0.085 


pVectorFFNs 


-0.085 




SigmaFFPd 


-0.23 


SigmaFFNd 


0.84 




SigmaFFPu 


0.84 


SigmaFFNu 


-0.23 


Tensor form factor 


SigmaFFPs 


-0.046 


SigmaFFNs 


-0.046 





References 



[1] Planck Collaboration, P. Ade et al, "Planck 2013 results. XVI. Cosmological 
parameters," 1303 . 5076. 

[2] R. Bernabei et al, "New results from DAMA/LIBRA," Eur. Phys. J. C67 (2010) 
39-49, 1002.1028. 

[3] CDMS Collaboration, R. Agnese et al, "Dark Matter Search Results Using the 
Silicon Detectors of CDMS 11," Phys. Rev. Lett. (2013) 1304.4279. 

[4] C. E. Aalseth et al, "Search for an Annual Modulation in a P-type Point Contact 
Germanium Dark Matter Detector," Phys. Rev. Lett. 107 (2011) 141301, 
1106.0650. 

[5] G. Angloher, M. Bauer, 1. Bavykina, A. Bento, C. Bucci, et al., "Results from 730 
kg days of the CRESST-11 Dark Matter Search," Eur.Phys.J. C72 (2012) 1971, 
1109.0702. 

[6] XENONIOO Collaboration, E. Aprile et al, "Dark Matter Results from 225 Live 
Days of XENONIOO Data," 1207.5988. 

[7] COUPP Collaboration, E. Behnke et al, "Improved Spin-Dependent WIMP 
Limits from a Bubble Chamber," Science 319 (2008) 933-936, 0804.2886. 

[8] S. Archambault et al, "Dark Matter Spin-Dependent Limits for WIMP 

Interactions on 19-F by PICASSO," Phys. Lett. B682 (2009) 185-192, 0907.0307. 

[9] Fermi-LAT Collaboration, M. Ackermann et al., "Constraining Dark Matter 
Models from a Combined Analysis of Milky Way Satellites with the Fermi Large 
Area Telescope," Phys. Rev. Lett. 107 (2011) 241302, 1108.3546. 

[10] C. Weniger, "A Tentative Gamma-Ray Line from Dark Matter Annihilation at the 
Fermi Large Area Telescope," JCAP 1208 (2012) 007, 1204.2797. 

[11] Fermi LAT Collaboration, M. Ackermann et al., "Measurement of separate 
cosmic-ray electron and positron spectra with the Fermi Large Area Telescope," 
Phys. Rev. Lett. 108 (2012) 011103, 1109.0521. 

[12] PAMELA Collaboration, O. Adriani et al., "An anomalous positron abundance 
in cosmic rays with energies 1.5-100 GeV," Nature 458 (2009) 607-609, 
0810.4995. 

[13] AMS Collaboration, M. Aguilar et al., "First Result from the Alpha Magnetic 
Spectrometer on the International Space Station: Precision Measurement of the 
Positron Fraction in Primary Cosmic Rays of 0.5-350 GeV," Phys. Rev. Lett. 110 
(2013), no. 14, 141102. 

[14] O. Adriani et al, "A new measurement of the antiproton-to-proton flux ratio up to 
100 GeV in the cosmic radiation," Phys. Rev. Lett. 102 (2009) 051101, 0810.4994. 



[15] Super-Kamiokande Collaboration, T. Tanaka et ai, "An Indirect Search for 
WIMPs in the Sun using 3109.6 days of upward-going muons in 
Super-Kamiokande," Astrophys.J. 742 (2011) 78, 1108.3384. 

[16] J. Zornoza, "Search for Dark Matter in the Sun with the ANTARES Neutrino 
Telescope in the CMSSM and mUED frameworks," 1204.5290. 

[17] IceCube Collaboration, R. Abbasi et ai, "Multi-year search for dark matter 

annihilations in the Sun with the AMANDA-II and IceCube detectors," Phys.Rev. 
D85 (2012) 042002, 1112.1840. 

[18] ATLAS Collaboration, G. Aad et ai, "Observation of a new particle in the search 
for the Standard Model Higgs boson with the ATLAS detector at the LHC," 
Phys.Lett. B716 (2012) 1-29, 1207.7214. 

[19] CMS Collaboration, S. Chatrchyan et ai, "Observation of a new boson at a mass 
of 125 GeV with the CMS experiment at the LHC," Phys.Lett. B716 (2012) 
30-61, 1207.7235. 

[20] ATLAS Collaboration, G. Aad et ai, "Search for dark matter candidates and 
large extra dimensions in events with a photon and missing transverse momentum 
in pp collision data at y/s = 7 TeV with the ATLAS detector," Phys.Rev. Lett. 110 
(2013) 011802, 1209.4625. 

[21] ATLAS Collaboration, G. Aad et ai, "Search for dark matter candidates and 
large extra dimensions in events with a jet and missing transverse momentum with 
the ATLAS detector," JHEP 1304 (2013) 075, 1210.4491. 

[22] CMS Collaboration Collaboration, S. Chatrchyan et ai, "Search for dark 
matter and large extra dimensions in mono jet events in pp collisions at ^/s = 7 
TeV," JHEP 1209 (2012) 094, 1206.5663. 

[23] P. Gondolo et ai, "DarkSUSY: Computing supersymmetric dark matter properties 
numerically," JCAP 0407 (2004) 008, astro-ph/0406204. 

[24] H. Baer, A. Belyaev, T. Krupovnickas, and J. O'Farrill, "Indirect, direct and 
coUider detection of neutrahno dark matter," JCAP 0408 (2004) 005, 
hep-ph/0405210. 

[25] A. Arbey and F. Mahmoudi, "Superlso Relic: A Program for calculating relic 
density and flavor physics observables in Supersymmetry," 
Comput.Phys.Commun. 181 (2010) 1277-1292, 0906.0369. 

[26] A. Semenov, "LanHEP - a package for the automatic generation of Feynman rules 
in field theory. Version 3.0," Comput. Phys. Commun. 180 (2009) 431-454, 
0805.0555. 

[27] R. Barbieri, L. J. Hall, and V. S. Rychkov, "Improved naturalness with a heavy 
Higgs: An Alternative road to LHC physics," Phys.Rev. D74 (2006) 015007, 
hep-ph/0603188. 



[28] G. Belanger, K. Kannike, A. Pukhov, and M. Raidal, "Impact of 

semi- annihilations on dark matter phenomenology - an example of Z^r symmetric 
scalar dark matter," JCAP 04 (2012) 010, 1202.2962. 

[29] U. Ellwanger and C. Hugonie, "NMHDECAY 2.0: An Updated program for 
sparticle masses, Higgs masses, couplings and decay widths in the NMSSM," 
Comput.Phys.Commun. 175 (2006) 290-303, hep-ph/0508022. 

[30] U. Ellwanger and C. Hugonie, "NMSPEC: A Fortran code for the sparticle and 
Higgs masses in the NMSSM with GUT scale boundary conditions," 
Comput.Phys.Commun. 177 (2007) 399-407, hep-ph/0612134. 

[31] A. Djouadi, J.-L. Kneur, and G. Moultaka, "SuSpect: A Fortran code for the 
supersymmetric and Higgs particle spectrum in the MSSM," Comput. Phys. 
Commun. 176 (2007) 426-455, hep-ph/0211331. 

[32] J. Lee, M. Carena, J. Ellis, A. Pilaftsis, and C. Wagner, "CPsuperH2.3: an 
Updated Tool for Phenomenology in the MSSM with Explicit CP Violation," 
Comput.Phys.Commun. 184 (2013) 1220-1233, 1208.2212. 

[33] G. Belanger, F. Boudjema, and A. Pukhov, micrOMECAs: A Package for 

Calculation of Dark Matter Properties in Ceneric Model of Particle Interaction, 
ch. 12, pp. 739-790. World Scientific, 2012. 

http : //www . worldscientif ic . com/doi/pdf /lO . 1142/9789814390163_0012. 

[34] G. Belanger, A. Pukhov, and G. Servant, "Dirac Neutrino Dark Matter," JCAP 
0801 (2008) 009, 0706.0526. 

[35] V. Barger, P. Langacker, M. McCaskey, M. Ramsey-Musolf, and G. Shaughnessy, 
"Complex Singlet Extension of the Standard Model," Phys. Rev. D79 (2009) 
015018, 0811.0393. 

[36] G. Belanger, K. Kannike, A. Pukhov, and M. Raidal, "Z3 Scalar Singlet Dark 
Matter," 1211.1014. 

[37] M. L. Graesser, I. M. Shoemaker, and L. Vecchi, "Asymmetric WIMP dark 
matter," J HEP 1110 (2011) 110, 1103.2771. 

[38] H. Iminniyaz, M. Drees, and X. Chen, "Relic Abundance of Asymmetric Dark 
Matter," JC^P 1107 (2011) 003, 1104.5548. 

[39] U. Ellwanger and P. Mitropoulos, "Upper Bounds on Asymmetric Dark Matter 
Self Annihilation Cross Sections," JCAP 1207 (2012) 024, 1205.0673. 

[40] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, "micrOMEGAs: Version 
1.3," Comput. Phys. Commun. 174 (2006) 577-604, hep-ph/0405253. 

[41] G. Belanger, B. Dumont, U. Ellwanger, J. Gunion, and S. Kraml, "Status of 
invisible Higgs decays," 1302 . 5694. 

[42] P. P. Giardino, K. Kannike, I. Masina, M. Raidal, and A. Strumia, "The universal 
Higgs fit," 1303.3570. 



[43] T. Hambye, "Hidden vector dark matter," JHEP 0901 (2009) 028, 0811.0172. 

[44] F. D'Eramo and J. Thaler, "Semi-annihilation of Dark Matter," JHEP 1006 
(2010) 109, 1003.5912. 

[45] M. Hindmarsh and O. Philipsen, "WIMP dark matter and the QCD equation of 
state," Phys.Rev. D71 (2005) 087302, hep-ph/0501232. 

[46] K. A. Olive, D. N. Schramm, and G. Steigman, "Limits on New Superweakly 
Interacting Particles from Primordial Nucleosynthesis," Nucl. Phys. B180 (1981) 
497. 

[47] M. Srednicki, R. Watkins, and K. A. Olive, "Calculations of relic densities in the 
early universe," Nucl. Phys. B310 (1988) 693. 

[48] G. Hinshaw, D. Larson, E. Komatsu, D. Spergel, C. Bennett, et al, "Nine- Year 
Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological 
Parameter Results," 1212.5226. 

[49] C. E. Yaguna, "Large contributions to dark matter annihilation from three-body 
final states," Phys.Rev. D81 (2010) 075024, 1003.2730. 

[50] A. Goudelis, B. Herrmann, and O. Stl, "Dark matter in the Inert Doublet Model 
after the discovery of a Higgs-like boson at the LHC," 1303.3010. 

[51] L. Lopez Honorez and C. E. Yaguna, "The inert doublet model of dark matter 
revisited," JHEP 09 (2010) 046, 1003.3125. 

[52] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, "micrOMEGAs2.0: A 
program to calculate the relic density of dark matter in a generic model," Comput. 
Phys. Commun. 176 (2007) 367-382, hep-ph/0607059. 

[53] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, "micrOMEGAs: A 

program for calculating the relic density in the MSSM," Comput. Phys. Commun. 
149 (2002) 103-120, hep-ph/01 12278. 

[54] IceCube Collaboration, R. Abbasi et ai, "Search for Neutrinos from Annihilating 
Dark Matter in the Direction of the Galactic Center with the 40-String IceCube 
Neutrino Observatory," 1210.3557. 

[55] M. Boliev, S. Demidov, S. Mikheyev, and O. Suvorova, "Search for muon signal 
from dark matter annihilations in the Sun with the Baksan Underground 
Scintillator Telescope for 24.12 years," 1301.1138. 

[56] IceCube Collaboration, C. H. Ha, "The First Year IceCube-DeepCore Results," 
J.Phys.Conf.Ser. 375 (2012) 052034, 1201.0801. 

[57] A. Gould, "Resonant Enhancements in WIMP Capture by the Earth," 
Astrophys.J 321 (1987) 571. 

[58] M. Asplund, N. Grevesse, and J. Sauval, "The Solar chemical composition," 
Nucl.Phys. A777 (2006) 1-4, astro-ph/0410214. 

An 



[59] W. McDonough, Compositional Model for the Earth's Core, vol. 2, p. 547568. 
Elsevier Ltd, 2003. 

[60] K. Griest and D. Seckel, "Cosmic Asymmetry, Neutrinos and the Sun," Nucl. 
Phys. B283 (1987) 681. 

[61] A. Gould, "WIMP distribution in and evaporation from the Sun," Astrophys.J. 
321 (1987) 560. 

[62] A. Gould, "Evaporation of WIMPs with arbitrary cross-sections," Astrophys. J. 
(1989). 

[63] M. Cirelli et ai, "Spectra of neutrinos from dark matter annihilations," Nucl. 
Phys. B727 (2005) 99-138, hep-ph/0506298. 

[64] G. Belanger, M. Kakizaki, E. K. Park, S. Kraml, and A. Pukhov, "Light mixed 
sneutrinos as thermal dark matter," JCAP 1011 (2010) 017, 1008.0580. 

[65] A. E. Erkoca, M. H. Reno, and I. Sarcevic, "Muon Fluxes From Dark Matter 
Annihilation," Phys. Rev. D80 (2009) 043514, 0906.4364. 

[66] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, "Dark matter direct 
detection rate in a generic model with micrOMEGAs2.1," Comput. Phys. 
Commun. 180 (2009) 747-767, 0803.2360. 

[67] G. Belanger et ai, "Indirect search for dark matter with micrOMEGAs2.4," 
Comput. Phys. Commun. 182 (2011) 842-856, 1004.1092. 

[68] F. Boudjema, A. Semenov, and D. Temes, "Self-annihilation of the neutralino dark 
matter into two photons or a Z and a photon in the MSSM," Phys. Rev. D72 
(2005) 055024, hep-ph/0507127. 

[69] N. Baro, F. Boudjema, and A. Semenov, "Automatised full one-loop 

renormalisation of the MSSM. I. The Higgs sector, the issue of tan(beta) and 
gauge invariance," Phys. Rev. D78 (2008) 115003, 0807.4668. 

[70] N. Baro and F. Boudjema, "Automatised full one-loop renormalisation of the 
MSSM II: The chargino-neutralino sector, the sfermion sector and some 
applications," Phys. Rev. D80 (2009) 076010, 0906.1665. 

[71] G. Chalons and A. Semenov, "Loop-induced photon spectral lines from neutralino 
annihilation in the NMSSM," JHEP 1112 (2011) 055, 1110.2064. 

[72] H. Zhao, "Analytical Models For Galactic Nuclei," Mon. Not. Roy. Astron. Soc. 
278 (1996) 488-496, astro-ph/9509122. 

[73] J. Lavalle, J. Pochon, P. Salati, and R. Taillet, "dumpiness of Dark Matter and 
Positron Annihilation Signal: Computing the odds of the Galactic Lottery," 
Astron. Astrophys. 462 (2007) 827-848, astro-ph/0603796. 

[74] A. Bottino, F. Donato, N. Fornengo, and S. Scopel, "Implications for relic 
neutralinos of the theoretical uncertainties in the neutralino nucleon 
cross-section," Astropart.Phys. 13 (2000) 215-225, hep-ph/9909228. 



41 



[75] J. R. Ellis, K. A. Olive, and C. Savage, "Hadronic Uncertainties in the Elastic 
Scattering of Supersymmetric Dark Matter," Phys.Rev. D77 (2008) 065026, 
0801.3656. 

[76] A. Thomas, P. Shanahan, and R. Young, "Strangeness in the nucleon: what have 
we learned?," Nuovo Cim. C035N04 (2012) 3-10, 1202.6407. 

[77] Particle Data Group Collaboration, J. Beringer et ai, "Review of Particle 
Physics (RPP)," Phys.Rev. D86 (2012) 010001. 

[78] JLQCD Collaboration, H. Ohki et ai, "Nucleon strange quark content from 
Nf = 2 + 1 lattice QCD with exact chiral symmetry," 1208.4185. 

[79] M. Engelhardt, "Strange quark contributions to nucleon mass and spin from 
lattice QCD," Phys.Rev. D86 (2012) 114510, 1210.0025. 

[80] MILC Collaboration, W. Freeman and D. Toussaint, "The intrinsic strangeness 
and charm of the nucleon using improved staggered fermions," 1204.3866. 

[81] R. Young and A. Thomas, "Octet baryon masses and sigma terms from an SU(3) 
chiral extrapolation," Phys.Rev. D81 (2010) 014503, 0901.3310. 

[82] S. Durr, Z. Fodor, T. Hemmert, C. Hoelbling, J. Prison, et ai, "Sigma term and 
strangeness content of octet baryons," Phys.Rev. D85 (2012) 014509, 1109.4265. 

[83] R. Horsley, Y. Nakamura, H. Perlt, D. Pleiter, P. Rakow, et ai, "Hyperon sigma 
terms for 2+1 quark flavours," Phys.Rev. D85 (2012) 034506, 1110.4971. 

[84] A. Semke and M. Lutz, "Strangeness in the baryon ground states," Phys.Lett. 
B717 (2012) 242-247, 1202.3556. 

[85] P. Shanahan, A. Thomas, and R. Young, "Sigma terms from an SU(3) chiral 
extrapolation," 1205 . 5365. 

[86] X.-L. Ren, L. Ceng, J. M. Camalich, J. Meng, and H. Toki, "Octet baryon masses 
in next-to-next-to-next-to-leading order covariant baryon chiral perturbation 
theory," 1209.3641. 

[87] P. Junnarkar and A. Walker-Loud, "The Scalar Strange Content of the Nucleon 
from Lattice QCD," 1301.1114. 

[88] G. Bali, P. Bruns, S. Collins, M. Deka, B. Glasle, et ai, "Nucleon mass and sigma 
term from lattice QCD with two light fermion flavors," Nucl.Phys. B866 (2013) 
1-25, 1206.7034. 

[89] A. Belyaev, N. D. Christensen, and A. Pukhov, "CalcHEP 3.4 for collider physics 
within and beyond the Standard Model," 1207 . 6082. 

[90] A. Djouadi, "The Anatomy of electro-weak symmetry breaking. I: The Higgs 
boson in the standard model," Phys.Rept. 457 (2008) 1-216, hep-ph/0503172. 



AO 



[91] A. Djouadi, "The Anatomy of electro-weak symmetry breaking. II. The Higgs 
bosons in the minimal supersymmetric model," Phys.Rept. 459 (2008) 1-241, 
hep-ph/0503173. 

[92] A. Djouadi, J. Kalinowski, and M. Spira, "HDECAY: A Program for Higgs boson 
decays in the standard model and its supersymmetric extension," 
Comput.Phys.Commun. 108 (1998) 56-74, hep-ph/9704448. 

[93] P. Baikov and K. Chetyrkin, "Higgs Decay into Hadrons to Order alpha**5(s)," 
Phys. Rev. Lett. 97 (2006) 061803, hep-ph/0604194. 

[94] P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein, and K. E. Williams, 

"HiggsBounds 2.0.0: Confronting Neutral and Charged Higgs Sector Predictions 
with Exclusion Bounds from LEP and the Tevatron," Comput.Phys.Commun. 182 
(2011) 2605-2631, 1102.1898. 

[95] G. B. et al, "Implementaion of new physics models in micrOMEGAs." in 
preparation. 

[96] G. Belanger, N. D. Christensen, A. Pukhov, and A. Semenov, "SLHAplus: a 
library for implementing extensions of the standard model," 
Comput.Phys.Commun. 182 (2011) 763-774, 1008.0181. 

[97] M. S. Carena, M. Olechowski, S. Pokorski, and C. Wagner, "Electroweak 
symmetry breaking and bottom - top Yukawa unification," Nucl.Phys. B426 
(1994) 269-300, hep-ph/9402253. 

[98] D. M. Pierce, J. A. Bagger, K. T. Matchev, and R.-j. Zhang, "Precision 

corrections in the minimal supersymmetric standard model," Nucl. Phys. B491 
(1997) 3-67, hep-ph/9606211. 

[99] F. Boudjema and A. Semenov, "Measurements of the SUSY Higgs selfcouplings 
and the reconstruction of the Higgs potential," Phys. Rev. D66 (2002) 095007, 
hep-ph/0201219. 

[100] M. S. Carena, M. Quiros, and C. Wagner, "Effective potential methods and the 
Higgs mass spectrum in the MSSM," Nucl.Phys. B461 (1996) 407-436, 
hep-ph/9508343. 

[101] B. AUanach, C. Balazs, G. Belanger, M. Bernhardt, F. Boudjema, et al., "SUSY 
Les Houches Accord 2," Comput.Phys.Commun. 180 (2009) 8-25, 0801.0045. 

[102] F. Mahmoudi, "Superlso v2.3: A Program for calculating flavor physics 

observables in SupersjTumetry," Comput.Phys.Commun. 180 (2009) 1579-1613, 
0808.3144. 

[103] FlaviaNet Working Group on Kaon Decays Collaboration, M. Antonelli et 
al., "Precision tests of the Standard Model with leptonic and semileptonic kaon 
decays," 0801.1817. 

[104] A. Akeroyd and F. Mahmoudi, "Constraints on charged Higgs bosons from 
D{s)^ /i±z/ and D{s)^ r±z/," JHEP 0904 (2009) 121, 0902.2393. 

A'i 



Heavy Flavor Averaging Group Collaboration, Y. Amhis et a/., "Averages of 
B-Hadron, C-Hadron, and tau-lepton properties as of early 2012," 1207. 1158. 

G. Belanger, F. Boudjema, C. Hugonie, A. Pukhov, and A. Semenov, "Relic 
density of dark matter in the NMSSM," JCAP 0509 (2005) 001, hep-ph/0505142. 

T. Halin and M. Perez- Victoria, "Automatized one loop calculations in 
four-dimensions and D- dimensions," Comput.Phys.Commun. 118 (1999) 153-165, 
hep-ph/9807565. 

F. Mahmoudi, S. Heinemeyer, A. Arbey, A. Bharucha, T. Goto, et ai, "Flavour 
Les Houches Accord: Interfacing Flavour related Codes," Comput.Phys.Commun. 
183 (2012) 285-298, 1008.0762. 

B. C. Allanach, "SOFTSUSY: a program for calculating supersymmetric spectra," 
Comput. Phys. Commun. 143 (2002) 305-331, hep-ph/0104145. 

W. Porod, "SPheno, a program for calculating supersymmetric spectra, SUSY 
particle decays and SUSY particle production at e+ e- colliders," Comput. Phys. 
Commun. 153 (2003) 275-315, hep-ph/0301101. 

F. E. Paige, S. D. Protopopescu, H. Baer, and X. Tata, "ISAJET 7.69: A Monte 
Carlo event generator for pp, anti-p p, and e+e- reactions," hep-ph/0312045. 

A. Djouadi, M. Muhlleitner, and M. Spira, "Decays of supersymmetric particles: 
The Program SUSY-HIT (SUspect-SdecaY-Hdecay-InTerface)," Acta Phys.Polon. 
B38 (2007) 635-644, hep-ph/0609292. 

G. Degrassi, P. Gambino, and P. Slavich, "SusyBSG: A Fortran code for 
BR{B X(s)7 in the MSSM with Minimal Flavor Violation," 
Comput.Phys.Commun. 179 (2008) 759-771, 0712.3265. 



AA 



